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A framework was recently introduced to generalize black hole mechanics by 
replacing stationary event horizons with isolated horizons. That framework 
is significantly extended. The extension is non-trivial in that not only do the 

OO ' 

boundary conditions now allow the horizon to be distorted and rotating, but 
■ also the subsequent analysis is based on several new ingredients. Specifically, 

although the overall strategy is closely related to that in the previous work, 
the dynamical variables, the action principle and the Hamiltonian framework 
*q \ are all quite different. More importantly, in the non-rotating case, the first law 

is shown to arise as a necessary and sufficient condition for the existence of 
^Jy a consistent Hamiltonian evolution. Somewhat surprisingly, this consistency 

condition in turn leads to new predictions even for static black holes. To 
• i-H . complement the previous work, the entire discussion is presented in terms of 

^ \ tetrads and associated (real) Lorentz connections. 



I. INTRODUCTION 



The zeroth and first laws of black hole mechanics refer to equilibrium situations and 
small departures therefrom. The standard treatments [|T]-|§] restrict themselves to stationary 
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space-times admitting event horizons and small perturbations from stationarity. While this 
simple idealization is a natural starting point, from physical considerations it seems overly 
restrictive. (See ||0] and especially [|J for a detailed, critical discussion.) A framework, 
which is tailored to more realistic physical situations was introduced in |J and the zeroth 
and first laws were extended to it in |7|-£|. This analysis generalizes black hole mechanics in 
two directions. First, the notion of event horizons is replaced by that of 'isolated horizons'. 
While the former can only be defined retroactively, requiring access to the entire space-time 
history, the latter can be defined quasi-locally. Second, the underlying space-time need not 
admit any Killing field; isolated horizons need not be Killing horizons. The static event 



horizons normally used in black hole mechanics |ffl-|,lC] and the cosmo logical horizons in de 



Sitter space-times [(TT] are all special cases of isolated horizons. Furthermore, since space- 
times can now admit gravitational and matter radiation, there is a large class of other 
examples. 

The framework developed in H for generalizing black hole mechanics was based on two 
restrictive assumptions. First, only undistorted, non- rotating horizons were considered. That 
is, the boundary conditions used in || implied that the intrinsic 2-metric of the horizon is 
spherically symmetric^ and that the imaginary part of the Weyl tensor component ^2 " 
which encodes the angular momentum — vanishes. Second, while a rather general class of 
matter fields was allowed, it was assumed that the only relevant charges — i.e., hair — are 
the standard electromagnetic ones. The second assumption was weakened in H9j,|l2]l which 
allowed dilaton couplings and Yang-Mills fields. In this paper, we allow for distortion and 
more general matter sources. Distortion plays an important role in several astrophysical 
situations, e.g., in problems involving black holes immersed in external fields, or surrounded 
by matter rings, and especially in the problem of black hole collisions. Post-Newtonian 
considerations suggest that, during black hole coalescence, individual horizons are distorted 
due to the Coulomb attraction even in the regime in which the black holes are sufficiently 
far from each other for the gravitational radiation falling into their horizons to be negligible. 
This phenomenon is also seen in numerical simulations. 

The extensions PJl2]| which incorporated dilatonic and Yang-Mills charges did not involve 



a significant generalization of the basic framework developed in ||. The present paper, on 
the other hand, does. We begin with substantially weaker boundary conditions (formulated 
in terms of (real) tensor fields rather than the spinors used in H), and show that they imply 
constancy of surface gravity (and electro-static potential) on the horizon. This property turns 
out to be necessary and sufficient for the usual action principle of tetrad gravity to continue 
to be valid in presence of isolated horizons. The action leads to a covariant phase space, 
constructed from solutions to the field equations. Ref ||, by contrast, used the canonical 
phase space based on spinorial variables which is tailored for quantization but which contains 



1 However, it allowed space-times, such as the Robinson- Trautman solutions, in which there is no 



space-time Killing field whatsoever in any neighborhood of the isolated horizon [15]. 
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technical complications that are unnecessary to the classical mechanics of isolated horizons. 
Up to this point, distortion and rotation are both incorporated. However, in the last step, i.e., 
in the discussion of the first law, we restrict ourselves to non-rotating horizons. (Rotation is 
incorporated in [0.) 

To formulate the first law one must first define the energy associated with any iso- 
lated horizon A. Since there can be radiation in the spacetime outside isolated horizons, 
the ADM energy Eadm is not a good measure of [0,[8| . Instead, as in the work of Brown 
and York |14| , the strategy is to define the energy of the horizon using a Hamiltonian frame- 
work. Experience with the phase space formulation of general relativity suggests that, in 
the presence of boundaries, the Hamiltonian H t generating time-translation along a suitable 
vector field t a acquires surface terms. The idea is to define as the surface term at A in 
the Hamiltonian. The key issue then is that of selecting the 'appropriate' time-translation 
t a . Since one expects the volume term in the expression of H t to be a linear combination 
of constraints and thus vanish when evaluated on solutions, the problem reduces to that of 
specifying the boundary values of t a (or, equivalently, of the lapse and shift fields). The 
conditions at infinity are obvious and, in any case, will not affect the surface term at A. 
Thus, we need to focus only on the boundary value of f on A. 

In the non- rotating case it is clear that, at the horizon, t a should be proportional to the 
null normal to A. However, our boundary conditions do not select the null normal uniquely; 
there is a freedom to rescalef] the normal by a constant (on A) which can vary from one 
space-time to another. This freedom is physically important because, amongst other things, 
the surface gravity is sensitive to it. Suppose we fix this freedom by tying the boundary 
value of t a to fields on the horizon, e.g., by demanding that surface gravity be a specific 
function of the area and charges. We can then ask whether the time-evolution generated by 
this t a preserves the symplectic structure. It turns out that the answer is not always in the 
affirmative. On the horizon, the evolution vector field t a and the electromagnetic potential 
A a have to be tied to the horizon parameters appropriately. These conditions impose a 
constraint on the surface gravity K(t) and the electric potential . Somewhat surprisingly, 
the constraint is precisely the first law SE^ = /87rG) 5aA + $(t) 5Qa- Thus the evolution 
defined by t a is Hamiltonian if and only if the first law holds.^ In this sense, the first law is 
even more fundamental than it is generally taken to be. Conceptually, this is perhaps the 



2 This is not surprising since this freedom exists already on Killing horizons. If the space-time 
is asymptotically flat and admits a static Killing field globally, one can eliminate this freedom 
by restricting oneself to that Killing field which is unit at infinity. However, this strategy is not 
available if there is radiation in the exterior region, or, as in the static solutions representing 
distorted black holes, the metric fails to be asymptotically flat. 

3 In the undistorted context, while this role of the first law was known to the authors of R§, its 



importance was not fully appreciated. The importance was noticed independently in [16] and used 



effectively in [03] to extract physical information on spherical black holes with Yang- Mills hair. 
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most striking feature of the present framework. 

The requirement that the first law hold is not sufficient to fix t a uniquely. Although 
every t a must be a null normal to A, the rescaling of t a from one spacetime to another can 
depend on the horizon parameters and the first law does not fully determine this parameter 
dependence. There is an infinite family of parameter-dependent vector fields t a each defining 
a consistent Hamiltonian evolution and a horizon energy E f A . By contrast, at infinity all 
these vector fields must tend to a time-translation of an universal flat metric, used in the 
construction of the phase space. Hence, there is only the familiar 3-parameter freedom in 
the definition of E t ADM , associated with the choice of a rest-frame. Furthermore, in any one 
space-time, we can eliminate it by simply going to the rest frame and thus extract the total 
mass Madm of the system. Although not necessary for mechanics, it is natural to ask if one 
can define a similar notion of mass of isolated horizons. The answer is in the affirmative 
in the Einstein-Maxwell theory. Let us require that t a should not only lead to a consistent 
Hamiltonian evolution but agree, on static solutions, with the static Killing field which is 
unit at infinity. Then the horizon value of t a is uniquely determined for all space-times 
in the phase-space. There is a preferred notion of time-translation, say t a Q . We can set 
Ma. — -^a an d regard Ma as the mass of an isolated horizon A. In the earlier work on 
non-distorted horizons |7|-p,[Ti2|] , the discussion of the first law was carried out only in the 
context of these preferred evolution vector fields t a D . That derivation is more closely tied 
to the traditional discussion of the laws in the static context but is not necessary from the 
more general perspective of isolated horizons. Nonetheless, the availability of a canonical 
definition of the mass Ma is useful for other applications of this framework, e.g. to numerical 
relativity. 

The paper is organized as follows. In Section |H], we specify the boundary conditions defin- 
ing general isolated horizons, allowing both distortion and rotation. We explain the role of 



these conditions, compare them with those used in [|q-|9|JT^| and work out their consequences, 
including the zeroth law. In Section [TIT] we introduce the Lagrangian framework based on 
tetrads and (real) Lorentz connections and in Section |TV| , the covariant phase space. The 
first law is discussed in Section [V|. Upto this point, the focus is on the Einstein-Maxwell the- 
ory (although incorporation of the dilaton is straightforward). In Section [VT] we extend the 
framework to incorporate Yang-Mills fields. The horizon mass is introduced in Section |VII 
and subtleties associated with the dilaton and Yang-Mills fields are discussed. For the conve- 
nience of readers who may not be familiar with distorted black holes, Appendix A presents 
a variety of examples and, for convenience of readers who work in the Newman-Penrose 
formalism, Appendix B summarizes the structure of isolated horizons in that framework. 

We have attempted to make this paper self-contained in terms of methodology and tech- 
nical details. However, the motivation behind isolated horizons and certain properties of our 
Hamiltonian are the same as those discussed in detail in ||. Since the inclusion of distortion 
does not add anything substantial to these issues, we have refrained from repeating that 
discussion in this paper. 



4 



II. STRUCTURE OF ISOLATED HORIZONS AND THE ZEROTH LAW 



In this section, we will introduce the basic definitions of isolated horizons and analyze 
their immediate consequences. The definitions will become progressively stronger. However, 
even the strongest boundary conditions are significantly weaker than requiring the horizon 
to be a Killing horizon for a local Killing vector field. By proceeding in steps, we will be 
able to keep track of the precise assumptions that are needed to obtain various results. Also, 
the availability of a hierarchy of definitions will be useful in other applications — such as 
numerical relativity and quantum gravity — which lie beyond the scope of the present paper. 

Let us begin by introducing some notation. Throughout this paper we assume that all 
manifolds and fields under consideration are smooth. Let M. be a 4-manifold equipped with 
a metric g a b of signature (—,+,+,+). Let A be a null hypersurface of (M.,g a b)- A future 
directed null normal to A will be denoted by £. (In this paper, the term 'null normal' will 
always refer to a future directed null normal.) Let q a b — g a b be the degenerate intrinsic 

metric on A.f\ A tensor q ab on A will be called an 'inverse' of q a b if it satisfies q ab q ac qbd — q c d- 
Thus q ab is unique only up to addition of terms of the form i^V^ for some vector field V 
tangential to A. The expansion 9^ of a specific null normal £ is defined by 9^ = q ab V ' a £b, 
where V a is the derivative operator compatible with g^. It is straightforward to check that 
is independent of the choice of q ab . With this structure at hand, we can now introduce 
our first definition. 

A. Non-expanding Horizons 

Definition 1: A 3-dimensional sub-manifold A of a space-time {Ai,g a b) is said to be a 
non-expanding horizon if it satisfies the following conditions: 

• (i) A is topologically S 2 x 1R and null; 

• (ii) The expansion 9^ of i vanishes on A for any null normal £; 

• (iii) All equations of motion hold at A and the stress-energy tensor T ab of matter fields 
at A is such that —T b a £ b is future directed and causal for any future directed null normal 
£. 

Note that if conditions (ii) and (iii) hold for one null normal £ they hold for all. 

The role of these conditions is as follows. The restriction on topology is geared to the 
structure of horizons that result from gravitational collapse. However, it can be weakened. 



4 Equalities which hold only at A will be denoted by '= ' and the pullback of a covariant index 
will be denoted by an arrow under that index, e.g. u a will denote the pullback of the 1-form uj a to 
A. 



5 



One can retain the requirement that the horizon have compact cross-sections but replace S 2 
by a manifold with higher genus. Our main analysis will extend to this case in a straight- 
forward manner. More generally, we can allow A to have non-compact cross-sections, as for 
example in the case of certain acceleration horizons. The results presented in this section, 
including our derivation of the zeroth law will go through. However, since such horizons 
extend to infinity, our Hamiltonian framework will have to be modified appropriately. Fi- 
nally, one could envisage incorporation of NUT charge. This extension would be even more 
subtle because, if all fields are smooth, A would be topologically S 3 and £ would provide 
a Hopf fibration. In this case, A would not admit any cross-sections which are everywhere 
transverse to I. This extension will be discussed elsewhere. 

Requirement (iii) is analogous to the dynamical conditions one imposes at infinity. While 
at infinity one requires that the metric (and other fields) approach a specific solution to the 
field equations (the 'classical vacuum') , at the horizon we only ask that the field equations 
be satisfied. The energy condition involved is very weak; it is implied by the (much stronger) 
dominant energy condition that is typically used. Thus, the first and the last conditions are 
quite tame. 

The key condition is (ii). It implies, in particular, that the horizon area is constant 
'in time' and thus incorporates the idea that the horizon is isolated without having to 
assume the existence of a Killing field. We will denote the area by a& and refer to 
defined by = as the horizon radius. All these conditions are satisfied on any 

Killing horizon (with 2-sphere sections) if gravity is coupled to physically reasonable matter 
(including perfect fluids, Klein-Gordon fields, Maxwell fields possibly with dilatonic coupling, 
Yang-Mills fields). 

Although the conditions in the definition are quite weak, they have surprisingly rich 
consequences. We will now discuss them in detail. In some of this analysis it is convenient 
to use a null-tetrad and the associated Newman-Penrose quantities (see Appendix B and 
references therein). Given a specific null normal field £ a to A, we can introduce a complex 
null vector field m a tangential to A and a real, future directed null field n a transverse to 
A so that the following relations hold: n ■ £ = — l,m • m = 1 and all other scalar products 
vanish. The quadruplet (£, n, m, m) constitutes a null-tetrad. There is of course an infinite 
number of null tetrads compatible with a given £, related to one another by restricted Lorentz 
rotations. Our conclusions will not be sensitive to this gauge-freedom. 

(a) Properties of £: Since £ a is a null normal to A, it is automatically twist free and geodesic. 
We will denote the acceleration of £ a by 

£ a V a £ b = Hl) £ h . (II.l) 

Note that the acceleration is a property not of the horizon A itself, but of a specific null 
normal to it: if we replace £ by £' = f£, then the acceleration changes via = fny) + C-zf ■ 
Since the twist of £ vanishes, the Raychaudhuri equation implies: 

Ctfw = K(q0 w - \e 2 {l) -aa- R ab t£ b 
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where a = m a m b V a £b is the shear of £ in the given null tetrad. Since 6^ vanishes on A, 
we conclude: o~ct + R a b£ a £ b = 0. The condition on the stress-energy tensor ensures that 
R a b£ a £ b = 8nGT a b£ a £ b is non-negative on A. Hence, we conclude: 

a = 0, and R ab t£ b = 0. (11.2) 

Thus, in particular, every null normal £ is free of expansion, twist and shear. 

(b) Conditions on the Ricci tensor: The second equation in (|II.2|) implies that the vector 
—R a b£ b is tangential to A. The energy condition and the field equations imply this vector 
must also be future causal. This means that R a b£ b must be proportional to £ a and hence, 
Rab^ b — 0. In the Newman-Penrose formalism this condition translates to: 

%o = lRab£ a £ b = and $oi = $io = ^R ab £ a m b = 0. (II.3) 

Since this statement is equivalent to R a b£ b = 0, it is gauge invariant, i.e. it does not depend 
upon the specific choice of null normal I and m. 

(c) A natural connection 1-form on A: Since £ is expansion, shear and twist free, there exists 
a one-form u a intrinsic to A such that 

VJ b = uj b (H.4) 

which in turn implies 

C £ q ab = 2VJb = . 

Thus, every null normal £ is a 'Killing field' of the degenerate metric on A. Furthermore, 
we will now show that 

2 e := im A m (H-5) 

is also invariantly defined. Since Ciq a b = 0, the space S of integral curves of £ is naturally 
equipped with a non-degenerate metric q ab (so that q a b on A is the pull-back of q ,). Denote 
by e afe the unique (up to orientation) unit alternating tensor on (S, q ). \ is the pull-back 
to A of e. Although £ is a 'Killing field' of the intrinsic horizon geometry, the space-time 
metric g a b need not admit a Killing field in any neighborhood of A. Robinson- Trautman 
metrics [15] and the Kastor-Traschen solutions [DJ provide explicit examples of this type. 



The 1-form uj will play an important role throughout this paper. It has an interesting 
geometrical interpretation. We can regard uj as a connection on the line bundle TA 1 - over 
A whose fibers are the 1-dimensional null normals to A. Under the rescalings £ i— > £ = f£, 
of the null normal £, it transforms via: 

io a ^ Cu a = uj a + V a In /. (II. 6) 
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(d) Induced Connection on A: Each metric submanifold M of M. admits a natural connection 
- one which is torsion-free and compatible with the induced metric on M. This connection 
is also canonically induced by the space-time connection V. However, since the induced 
metric q ab on A is degenerate, there exist infinitely many connections compatible with it. A 
general null sub-manifold inherits a unique (torsion- free) derivative operator T> from V if and 
only if its null normal £ satisfies y a £ b = 0. Therefore, the conditions imposed in Definition 
1 guarantee that every non-expanding horizon has a unique intrinsic derivative operator D. 
The action of P on a vector field X a tangent to A and on a 1-form r] a intrinsic to A is given 
by 

V a X b = VaX b and V aVb = VJjb ■ 

where X b and fj b are arbitrary extensions of X b and r\ b to the full space-time manifold Ai. 
It is easy to show that T> is independent of the extensions. 

The 1-form uj captures only part of the information in T>. The full connection Don A 



plays an important role in extracting physics in the strong field regime near A [IB]. However, 
it is not essential to the discussion of isolated horizon mechanics. 

(e) Conditions on the Weyl tensor: Let us begin with the definition of the Riemann tensor, 
[V a V& — V&V a ]X c = —2R abd c X d . If we set X c = £ c and pull back the indices a and 6, then 
using ( |I1.4j ), we obtain: 

[V a u b - V b u a ]t = -2R abd c £ d = -2C abd c £ d (11.7) 

where C abc d is the Weyl tensor. The last equality follows from R ab £ b = 0. Thus, if v is any 
1-form on A satisfying v ■ £ = 0, contracting the previous equation with v c we get 

C a M C V c £ d = 0. 

Let us choose a null tetrad and set v to be m or m. Then 

*o := C abcd £ a m b £ c m d = and *j := C abcd £ a m b £ c n d = C abcd £ a m b m c m d = , (II.8) 

where we have used the trace-free property of the Weyl tensor in the second equation. It is 
also clear that equations ( |11.8j ) are independent of which null normal £, and vector fields m 
and m we choose to construct the null tetrad; equation (|II.8|) is gauge invariant. 



(f) Curvature of uj: Let us contract ( |11.7| ) with n c and use £ a n a = —1. Then we have: 

2V [a uj b] = C abd c £ d n c = C ab£d £ c n d . (II.9) 

Expanding the Weyl tensor in terms of the ^'s, one obtains 

C abcd £ c n d = 4(Re [^ 2 })n [a l b] + 2^ 3 £ [a m b] + 2^ 3 £ [a m b] 

-2 : * 1 n [a m b] - 2^ in[a m b] + 4z(Im [^ 2 \)m [a m b] . (11.10) 
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where 



^2 := C abcd £ a m b m c n d and V 3 := C abcd £ a n b m c n d . (11.11) 

Substituting this expression into ( |II.9| ), pulling back on the two free indices and taking into 
account ( [11. 8[ ) and ( |11.5| ), we obtain 



du = 2(Im[^ 2 ]) 2 e. (11.12) 

This relation will play an important role in what follows. Note that, because \l/o an d V&i 
vanish on A, ^ 2 is gauge invariant. 

Remark: It is interesting to compare the structure of A with that of null-infinity X (in 
the usual conformal gauge, in which the conformal factor is chosen such that the null- 
normal to X is divergence-free). Both are null surfaces and can be regarded as 'line bundles' 
over a base space S of the integral curves of null normals. (For brevity, we will ignore a 
caveat concerning completeness of fibers.) The null normals are Killing fields of the intrinsic 
degenerate metric so that this metric is the pull-back to the 3-surface of a positive-definite 
metric on S. In both cases, the space-time connection induces an intrinsic connection on the 



3-surface ||19|| . These connections capture physically important information. However, there 
are a number of differences as well. Since X is constructed by a conformal completion, the 
conformal freedom permeates all geometric structures at X. In particular, given a physical 
space-time, the intrinsic metric and the derivative operator are known only up to conformal 
transformations. On the other hand, since X is at infinity, in some ways its structure is both 
more rigid and simpler. First, without loss of generality, we can assume that the metric 
on S is a 2-sphere metric; the issue of distortion is physically irrelevant at X. Second, the 
Weyl tensor vanishes identically at X and the curvature of the intrinsic connection captures 
non-trivial information about the next order space-time curvature. By contrast, at A only 
four components of the Weyl curvature vanish and four other components are coded in the 
curvature of the intrinsic connection on A. In spite of these differences, one can carry over 
some techniques from null infinity to extract physical information about isolated horizons. In 
particular, using the analogs of techniques which have been successful at X, one can introduce 



preferred cross sections of and Bondi-type expansions near A ]T8| 



This concludes our analysis of the consequences of the boundary conditions defining non- 
expanding horizons. Note that, even though i is a Killing field for the intrinsic, degenerate 
metric qw, on A, it is not an infinitesimal symmetry for other geometrical fields such as the 
intrinsic connection T> or components of the curvature tensor. In the next sub-section, we 
will make the structure more rigid by suitably restricting the choice of I. 



B. Weakly Isolated horizons 

The time-independence of the intrinsic metric q a b captures the idea that A is isolated in 
a suitable sense. While this condition has rich consequences, the resulting structure is still 
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not sufficient for physical applications. In particular, since £ can be rescaled by an arbitrary 
positive definite function, the acceleration Km is not necessarily constant on A. Therefore, 
we need to impose additional restrictions on the physical fields at A to establish the zeroth 
law. Since £ is already a symmetry of the intrinsic metric, it is natural to require it also be 
a symmetry of the 'extrinsic curvature'. However, the standard definition of the extrinsic 
curvature is not applicable to null surfaces. Nonetheless, given a null normal £, we can 
construct a tensor field K b := T> a £ b , defined intrinsically on A, which can be thought of as 
the analogue of the extrinsic curvature.^ Indeed, on a metric sub-manifold, if we replace £ 
by the unit normal, K b is precisely the extrinsic curvature. It is then natural to demand 
that, on an isolated horizon, K a h also be time-independent: CeK a b = 0. As a consequence 
of QH.4|) , this is equivalent to imposing Ceu = 0. 



Let us examine the above condition. As we will show at the end of this section, given 
a non-expanding horizon we can always find a null normal £ a which satisfies Lgu = 0. The 
behavior of this condition under rescalings of £ is complicated by the fact that uo itself depends 
upon the choice of null normal (see equation (|II.6| )). However, under a constant rescaling 
£ i — > £ — c£^ the connection 1-form u is unchanged. Therefore, if £ satisfies the condition 
CiU = 0, so does any £ related to £ by constant rescaling. This suggests we introduce an 
equivalence relation: Two future-directed null normals £ and £ belong to the same equivalence 
class [£] if and only if £ = c£ for some positive constant c. 

The above considerations lead us to the following definition: 

Definition 2: A weakly isolated horizon (A, [£]) consists of a non-expanding horizon A, 
equipped with an equivalence class [£] of null normals to it satisfying 

C e u = for all £&[£]. (11.13) 

As pointed out above, if this last equation holds for one £, it holds for all £ in [£]. 

A Killing horizon (with 2-sphere cross-sections) is automatically a weakly isolated hori- 
zon, (provided the matter fields satisfy the energy condition of Definition 1). Given a non- 
expanding horizon A, one can always find an equivalence class [£} of null- normals such that 



(A, [£]) is a weakly isolated horizon. However, condition ( |II.13| ) does not by itself single out 
the appropriate equivalence class [£}. As indicated in Section [11 D] , one can further strengthen 
the boundary conditions and provide a specific prescription to select the equivalence class 
[£} uniquely. However, for mechanics of isolated horizons, these extra steps are unnecessary. 



5 We are grateful to Thibault Damour for pointing out that K a b is called the Weingarten map and 
is analogous to extrinsic curvature. This comment suggested the above motivation for our condition 
on the connection 1-form u. For an alternate, and in a sense weaker, condition see the remark at 
the end of Section [ID. From the viewpoint of intrinsic structures on A discussed in Section II A, it 
is perhaps more natural to ask that £ be a symmetry of the full intrinsic connection T> (see section 
II Dj and ]l8fl -) However, this stronger condition is not necessary for the laws of mechanics discussed 



here. 
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In particular, our analysis will not depend on how the equivalence class [£] is chosen. The 
adverb 'weakly' in Definition 2 emphasizes this point. 

The condition ( |11.13| ) has several consequences which are relevant for this paper. 



(a) Surface gravity: In the case of Killing horizons Ak, surface gravity is defined as the 
acceleration of the Killing field £ normal to Ak- However, if Ak is a Killing horizon for £, 
it is also a Killing horizon for c£ for any positive constant c. Hence, surface gravity is not 
an intrinsic property of Ak, but depends also on the choice of a specific Killing field £. (Of 
course the result that the surface gravity is constant on Ak is insensitive to this rescaling 
freedom.) In asymptotically flat space-times admitting global Killing fields, this ambiguity is 
generally resolved by selecting a preferred normalization in terms of the structure at infinity. 
For example, in the static case, one requires the Killing field £ to be unit at infinity. However, 
in absence of a global Killing field or asymptotic flatness, this strategy does not work and 
we simply have to accept the constant rescaling freedom in the definition of surface gravity. 
In the context of isolated horizons, then, it is natural to keep this freedom. 

A weakly isolated horizon is similarly equipped with a preferred family [£] of null normals, 
unique up to constant rescalings. Therefore, it is natural to interpret Km as the surface 
gravity associated with I. Under the permissible rescalings £ i— > £ = c£, the surface gravity 
transforms via: k,^ = CKyy Thus, while u is insensitive to the rescaling freedom in [£], Km. 
captures this freedom fully. One can, if necessary, select a specific £ in [£] by demanding that 
Km be a specific function of the horizon parameters which are insensitive to this freedom, 
e.g., by setting Km = 1/2R&, where is the horizon radius (related to the horizon area 
a A via «A = 47ri? A ). 

(b) Zeroth law: The boundary conditions of Definition 2 allow us to define surface gravity Km 
of a weakly isolated horizon (A, [£]). We will now show that the surface gravity is constant 
on A. In other words, the zeroth law holds for weakly isolated horizons. 



Recall from ( |II.12|) that on a non-expanding horizon, aluj = 2 Im [^ 2 ] e f° r an y choice of 
null normal £. Since \ is the pull-back to A of the alternating tensor e on the space S (of 
orbits of £), clearly £ ■ \ = 0. Therefore, 

£ ■ duj = 

for every null normal £. In particular, on a weakly isolated horizon this equation holds for 
any £ e [£]. Moreover, each of these restricted null normals also satisfies 

= £ e u = d(£-u)+£- du 

Hence, we conclude: 

d{£ ■ lo) = d(Km) — , 



where we have used the definition ( [II. 1| ) of Km. Thus, surface gravity is constant on A. 



Although this proof of the zeroth law appears extremely simple, the argument is not as 
trivial as it might first appear since we have used a number of consequences of the boundary 
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conditions derived in section [LI Aj . In contrast to earlier derivations @JT3] we do not require 



the presence of a Killing field even in a neighborhood of A. Therefore the proof applies also 
to space-times such as the Robinson- Trautman solutions [[5] which do not admit a Killing 



field. Also, A need not be 'complete' — it may be of finite affine length with respect to 
any £ — and may not admit the analog of a 'bifurcate surface' on which the Killing field 
vanishes. Finally, the field equations are used rather weakly; we only need to assume that 
— (R a b — \Ro~ a b)£ b is a future directed causal vector. 

Surface gravity does not have a definite value on a weakly isolated horizon. The value of 
Km depends upon the choice of null normal £ G [£]. Since all the normals £ to A are future 
directed, the rescaling constant c is necessarily positive. Therefore, if the surface gravity is 
non-zero (respectively, zero) with respect to one £, it is non-zero (respectively, zero) with 
respect to any other £ G [£]. This rescaling freedom is the same as the one discussed above 
in the context of Killing horizons. 

We will conclude this sub-section with three remarks. 

i) Freedom in the choice of [£]: Given a non-expanding horizon A, it is natural to ask if 
one can always select an equivalence class [£] of null normals such that (A, [£]) is a weakly 
isolated horizon. As indicated earlier in this section, the answer is in the affirmative and, 
furthermore, there is a considerable freedom in the choice of [£]. Let us examine this issue 
in some detail. 

Since I ■ du = for any null normal £ to a non-expanding horizon, it follows that a null 
normal I satisfies L(io = if and only if dn^ = 0. Thus, to find a family [£] required in the 
definition of weakly isolated horizons, it is necessary and sufficient to find a null normal £ 
for which the surface gravity is constant. On a non-expanding horizon, the surface gravity 
transforms as follows. 

If I £ = f£, then k^ = fnm, + C e f. 

Hence, starting with any £, we can simply solve for / by requiring that k,^ be constant on 
A. The solution is not unique. If K£ is constant, given any non-zero function g satisfying 
Ltg = and a constant k, let us set 

/ = ge~ K W v + — 
K {e) 

where v satisfies C^v = 1. Then, we obtain an £ G" [£} for which = k. This is the 
only freedom if both and to be constant. Thus, each non-expanding horizon 

gives rise to an infinite family of weakly isolated horizons. Put differently, although one 
can easily obtain weakly isolated horizons from non-expanding ones by choosing appropriate 
null normals [£], a specific weakly isolated horizon carries much more information than the 
non-expanding horizon it comes from. At the end of Section [H D| , we will indicate how one 
can further strengthen the boundary conditions to give a prescription for selecting a specific 
[£]. However, the analysis of this paper does not depend on how this selection is made. 

ii) How does Definition 2 compare with that used in the undistorted, non-rotating case? 
As one would expect, the definition given in is significantly stronger. Furthermore, it 
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was tied to a foliation from the beginning. More precisely, it assumed that there exists a 
foliation to which u is normal, with uj = —ku, and it required that the expansion Re \p] of the 
null normal n to the leaves of the foliation be constant on A (see appendix [FJ for definitions 
of the NP spin coefficients). Although it was shown that the foliation is unique if it exists, 
the heavy reliance on the foliation right from the beginning made that definition less elegant 
and its invariant content less transparent. Also, since we now allow du to be non-zero and 
impose no conditions on Re \p], we can now incorporate rotation and distortion. 

iii) Alternate boundary conditions: In the definition of isolated horizons, we required 
L(io = 0, which in particular implies C^Km = Ci(l ■ u) = 0. Thus, the definition itself 
assured us that Km is 'time independent' and to prove the zeroth law we had to show that 
it is also independent of 'angles'. Could we have used another definition in which the 'time 
dependence' of Km was not explicitly required but followed from other conditions? The 
answer is in the affirmative: In place of Ctu = 0, we could have required that A admit 
a foliation on which the expansion Re [//] of n and the Newman- Penrose spin coefficient 7r 



which carries the angular momentum information are 'time independent' |R|. Then the 
field equations would have implied that Km, is time-independent. Furthermore, all results 
of this paper go through (and were in fact first obtained) with these modified boundary 
conditions. Note however that the new condition is neither weaker nor stronger than the 
one we used. Both require that tt be time independent. In addition, the present definition 
of isolated horizons requires that Km be time independent while the alternative definition 
would have required, instead, that Re \p] be time independent. In this paper, we chose the 
present definition because it can be stated without reference to a foliation. 



C. Electromagnetic Field 

We shall now describe the form of the electromagnetic field at an isolated horizon and 
introduce a partial gauge fixing at the horizon which will allow us to introduce the notion 
of an electric potential. In the next three sections — where we discuss the action, phase 
space and first law — we will assume that the only matter field present at the isolated 
horizon are Maxwell fields. However, as our discussion will make it clear, this restriction is 
made primarily for simplicity. The overall framework is rather general and can accommodate 
matter for which there exists a well-defined action principle and a (covariant) Hamiltonian 



framework. In particular, in Section [VT|, we describe how to extend the formalism to include 
Yang Mills Fields. 

The isolated horizon boundary conditions restrict matter primarily through conditions 
on the stress-energy tensor T^. Let us begin with T a ^l a l h = 0, a direct consequence of 
the boundary conditions and Raychaudhuri equation. (This restriction arises due to the 
fact that A is a non-expanding horizon; the subsequent stronger boundary conditions do 
not further constrain F). Although this condition is weak, it turns out to have interesting 
consequences on the form of the electromagnetic field, F, at A. The stress-energy tensor for 
electromagnetism is given in terms of the field strength F as 
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T ab = ^[F ac F 6 c - -a ab F cd F cd }. (11.14) 

47T 4 

Let us contract this expression with £ a £ b and examine the consequences for F. This gives 

= T ab £ a £ b =\£ a m b F ab \ 2 , (11.15) 

where, to obtain the last expression, we have used the anti-symmetry of F and the fact that 
the metric at the horizon can be expressed in terms of a null-tetrad as g ab = —2£^ a n b ^ + 
2m( a m b y An immediate consequence of QII.15D is that 



ab 



0. (11.16) 



In order to obtain a similar expression for *F recall that the stress energy tensor can be 
rewritten as T ab = — ^ r [*F ac *F fe c — ^g ab *F cd *F cd ] . Applying the same argument which led to 
( |11.16| ), we obtain a similar restriction on *F, namely 



Fab — 0. (11.17) 

These two restrictions tell us there is no flux of electromagnetic radiation across the horizon. 

It is straightforward to show that ( |11.16| ), ( |11.17| ) and the form of the metric at A place 
further restrictions on T ab : 

T ab £ a m b = T ab £ a m b = 
T ab m a m b = T ab m a m b = 0. 

The first two equations contain no new information since we already knew from general 
arguments (see equation ( |11.3| )), that $ 10 = and $ i — 0. However, the last two equations 
do place further restrictions on the stress energy tensor. Since the equations of motion are 
enforced at the boundary, we see immediately that they are equivalent to further restricting 
the Ricci tensor at the horizon by requiring: 

$02 := l -Rabm a m b = and $ 20 := ^R ab rn a m b = . (11.18) 

This result need not hold for general matter fields; it relies on the properties of the electro- 
magnetic stress-energy tensor ( 11.14 ). 

The next task is to define the electric and magnetic charges of the horizon. Since the 
horizon is an inner boundary of spacetime, the normal to a 2-sphere cross section of the 
horizon will naturally be inward pointing. Bearing this in mind, we define the electric and 
magnetic charges of the horizon as 

Qa :=-i-/ *F and P A := -— / F. (11.19) 

47T JSa 47T JSa 

For these definitions to be meaningful, the values of Qa and Pa should be independent of 
the cross section of the horizon We will now show that A being a non-expanding horizon 
guarantees this is the case. Let us first evaluate 
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C e F = £-dF + 



The first term on the right hand side vanishes due to Maxwell's equations on A, while the 
second term is zero due to the previous restriction on F, ( |II.16|) . Therefore we conclude 



that F is Lie dragged by t. An identical argument for *F leads to the analogous conclusion. 
Therefore we obtain 

£^F = and £ € *F = 0. (11.20) 



This result, along with QH.16|) and (|II.17|) , guarantees that Qa and Pa are independent of 



the choice of cross section Sa of the horizon. Note that this result was obtained using only 
the boundary conditions; equations of motion in the bulk are not needed. 

Finally, let us examine the remaining freedom in the electromagnetic field. The boundary 
conditions do not restrict F a i,n a m b and *F a i,n a m b at all. These components describe the 
electromagnetic radiation flowing along the horizon. Therefore, isolated horizon boundary 
conditions allow electromagnetic radiation arbitrarily close to — and even at — the horizon, 
provided none crosses it. 

So far we have confined ourselves to the field strengths F and *F. However, in the action 
principle and the Hamiltonian framework we have to consider also the Maxwell potential A. 
Now, if the magnetic charge is non-zero, either one has to allow 'wire singularities' in the 
vector potentials or regard A as a connection on a non-trivial [/(l)-bundle. (If we regard it 
as a connection on a H + -bundle, the magnetic charge is necessarily zero.) Since we wish to 
deal only with smooth fields, we will not allow 'wire singularities' in the potentials. If we 
work with bundles, the magnetic charge is quantized whence the space of histories has several 
disconnected components. Thus, in the first law, we will not be able to consider variations 
5 of fields with 5P 7^ 0. As far as mechanics of isolated horizons is concerned, there is 
essentially no loss of generality if we restrict ourselves to the case Pa — 0. Therefore, in the 
next three sections, while working with Maxwell fields, we will do so. As usual, our final 
results can be formally extended to the case of non- vanishing magnetic charge by performing 



a duality rotation on F. (As discussed in Section [VI], the situation is rather different in the 
case of Yang-Mills fields.) 

Recall that the first law in the Einstein-Maxwell case involves the electro-static potential 
$. In static space-times, one typically sets $ = — £ a A a where £ is the static Killing field and 
the gauge is chosen such that the vector potential A tends to zero at infinity and satisfies 
C^A = everywhere in space-time. We now need a strategy to define the electric potential 
$ without reference to a Killing field or infinity. To this end, we introduce the following 
definition: 

Definition 3: The electromagnetic potential A will be said to be in a gauge adapted to the 
weakly isolated horizon (A, [£]) if it satisfies 

C e A = 0. (11.21) 

Mathematically, this restriction is analogous to this condition £ = imposed on the 
gravitational field in Definition 2. However, while the condition on u is a physical restriction 
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on the form of the gravitational field at A, the condition on A is a gauge choice; it can 
always be imposed without physically constraining the electromagnetic field strength. 

Given an electromagnetic potential A in a gauge adapted to (A, [£]), we can now define 
the scalar potential at the horizon in an obvious fashion: 

$ w := -£ ■ A. 

The key question now is whether our boundary conditions are strong enough to ensure that 
is constant on A. Only then can we hope to use this notion of the scalar potential in the 
first law. Note that this question is rather similar to the one we asked of surface gravity Km 
in Section II B| . By using arguments completely analogous to those that led us to the zeroth 



law, we will now show that the answer to the present question is also in the affirmative. In 
a gauge adapted to the horizon, 

= C e A = fvF - d$ (/) . 

As we saw above, the boundary conditions imply £■¥ = (equation ( |II.16|) ). Hence, it 
follows immediately that is constant on the horizon. We can regard this result as the 
'electromagnetic part' of the zeroth law of isolated horizon mechanics. 

We will see in the next section that these zeroth laws play a key role in making the 
gravitational and the electromagnetic action principles well-defined in presence of isolated 
horizons. As with surface gravity Km, the functional dependence of $m on the horizon 
parameters varies with the choice of i G [£]. We will see in Section [V] that the Hamiltonian 
framework constrains these dependencies in an interesting fashion. 



D. Other definitions and remarks 

In this sub-section, we introduce two new definitions which are important to the general 
framework of isolated horizons. 

The first is concerned with rotation. From one's experience with the Newman- Penrose 
framework, one expects the gravitational contribution to angular momentum to be coded in 



the imaginary part of ^2- This expectation will be shown to be correct in ||13|| . Therefore, 
in the Einstein-Maxwell theory, we introduce the following definition: 

Definition 4: A weakly isolated horizon (A, [£}) will be said to be non-rotating if Im [^2] 
vanishes on A. 



If (Ak, [£]) is a Killing horizon and £ is a hypersurface orthogonal, time-like vector field near 
Ak, on physical grounds one would expect the horizon to be non- rotating. Is this expectation 
compatible with our definition? The answer is in the affirmative. For, in this case, one can 
show that B a i, := C ac t,d!; c !; d vanishes in the region where £ a is time-like. Hence, by continuity, 
it also vanishes on A forcing Im [^2] to vanish there. Similarly, if the space-time admits a 
hypersurface orthogonal, rotational Killing field (f a in a neighborhood of A Im [$2] again 
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vanishes on A. The definition is again compatible with one's intuition that the horizon 
should be non-rotating in this case. In this paper, while we allow for presence of rotation in 
the first four sections, we will restrict ourselves to non-rotating horizons in the proof of the 
first law in Section [V[ 

Finally, for completeness, let us introduce a stronger notion of 'isolation' by strengthening 
the boundary conditions of Definition 2. 

Definition 5: A weakly isolated horizon (A, [£}) is said to be isolated if 

[£ £l V}V = (11.22) 

for all vector fields V tangential to A and all £ G [£] 

As before, if any one £ satisfies this condition, so do all £ G [£}. However, unlike (|II.13|) , 



condition ( |11.22| ) is a genuine restriction in the sense that it can not always be met by a 



judicious choice of null normals. Generically it does suffice to single out the equivalence class 
[£] uniquely [0. In particular, in the Kerr family, the only [£] which satisfies ( 11.22 ) is the 



one containing constant multiples of the globally defined Killing field which is orthogonal to 
the horizon. Every Killing horizon is of course an isolated horizon. Thus, even though ( |11.22j ) 
is a stronger condition than ( |II.13|) , it is still very weak compared to conditions normally 
imposed. For most physical applications, e.g., to numerical relativity, it is appropriate to 
work with isolated horizons. For mechanics of isolated horizons, however, we can — and 
will — work with the larger class of weakly isolated horizons. 

The following consequences of the boundary conditions defining weakly isolated horizons, 
derived in this section, will play an important role in the subsequent discussion: 

a) Va£b = 0Ja£ b ] 

b) Km = £ a uJ a , the surface gravity defined by the null normal £ a , is constant on A; 

c) There is a natural (area) 2-form \ on A satisfying £ % \ = and 2 e a b£ b = 0. 

d) The electromagnetic potential A is chosen to satisfy £ ?A = and in this gauge the 

scalar potential $(£) := — £ a A a is constant on A; and 

e) The electromagnetic field satisfies £ a F ab = 0; and £ a *F ab = 0. 



III. ACTION 

In this paper we use the first order formulation of general relativity in terms of tetrads 
and connections. Since tetrads are essential to incorporate spinorial matter, it is natural to 
base the framework on tetrads from the beginning. The use of a first order formalism, on the 
other hand, is motivated primarily by mathematical simplicity. In the first order framework, 
the action and the Hamiltonians can be expressed entirely in terms of differential forms 
which significantly simplify the variational calculations. The previous paper || which dealt 
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with undistorted horizons used spinors and self-dual connections, while here we choose to 
use orthonormal tetrads and real, Lorentz connections. For analyzing mechanics of isolated 
horizons, there are two advantages to this. First, the Hamiltonian and symplectic structure 
are now manifestly real which simplifies evaluation of the boundary terms at the horizon. 
Second, the analysis can now be extended to other space-time dimensions in a straightforward 
manner. However, these simplifications, come with a price. Since, at present, the self dual 
variables appear to be indispensable for non-perturbative quantization, the results obtained 
here will have to be re-expressed in terms of self-dual variables in order to extend the analysis 
pOfl of the quantum horizon geometry and black hole entropy to include rotation. 



A. Preliminaries 

Let us begin with the first order action for Einstein-Maxwell theory on a 4-dimensional 
manifold M. which is topologically M x IR, where M is an oriented Riemannian 3-manifold 
without boundary (the complement of a compact set of) which is diffeomorphic to (the 
complement of a compact set of) IR 3 . Thus, topological complications of M, if any, are 
confined to a compact set. In this subsection we shall only give the relevant formulae. For 
details, see e.g. [^l]]. Our basic fields will consist of a triplet (e T a ,A aI J , A a ) defined on M. 



where denotes a co-tetrad, A a j J the gravitational (Lorentz) connection and A a the electro- 
magnetic connection. Here, lower case latin letters refer to the tangent space of A4 while 
the upper case letters /, J etc. refer to an internal four dimensional vector space V with a 
fixed metric rju of signature ( — h ++). The co-tetrad is an isomorphism between the 
tangent space T P (J\4) at any point p and the internal space V. Using it, we define a metric 
on M. by g a b := e^rjjj which also has signature ( — h ++). The Lorentz connection A a j J 
acts only on internal indices and defines a derivative operator 

D a ki := d a kj + A aI J kj , 

where d is a fiducial derivative operator which, as usual, will be chosen to be flat and torsion 
free. Finally, A a is the U(l) electromagnetic connection 1-form on M.. (As noted in Section 
II C| , we will assume that the magnetic charge is zero.) All fields will be assumed to be 
smooth and satisfy the standard asymptotic conditions at infinity. 
The 2-forms S /J 

S/j := ^e IJKL e K A e L 

constructed from the co-tetrads will play an important role throughout our calculations. In 
particular, the action for an asymptotically flat space-time (with no internal boundary) is 



given by (see e.g. [21 



S(e,A,A) = — / Z IJ AF IJ + —!— f J: IJ AA U -—I F A *F . (III. 1 ; 



16nG Jm 16nG j 

Too 
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Here F and F are the curvatures of the gravitational and electromagnetic connections A and 
A respectively: 

F/ = dAf J + Aj K A A K J and F = dA , 

*F ab = \e ab cd ¥ cd is the dual of F defined using ej, and is the time-like cylinder at 
infinity The boundary term at ensures the differentiability of the action. 

Let us briefly examine the equations of motion arising from the action. Varying the 
action with respect to the connection, one obtains 

L>£ = 0. 

This condition implies that the connection D defined by A has the same action on internal 
indices as the unique connection V compatible with the co-tetrad, i.e., satisfying V ' a e[ = 0. 
When this equation of motion is satisfied, the curvature F is related to the Riemann curvature 
R of V by 

771 IJ _ p cd I J 

f a b — K a b e c e d . 

Varying the action with respect to e a and taking into account the above relation between 
curvatures, one obtains Einstein's equations 

G ah = SnGTab = 2G (F ac F bd g cd - ^g ah F cd F cd ^) 

where G ab is the Einstein tensor and T ab the electromagnetic stress energy tensor. Finally, 
variation with respect to the electromagnetic connection, A, yields Maxwell's equation 

d*F = . 

B. Internal boundary A 

Let us now consider the variational principle for asymptotically flat histories which admit 
a weakly isolated horizon A as their internal boundary. The manifold Ai under consideration 
has an internal boundary A, topologically S 2 x 1R. As before, M. is topologically M x IR, 
where M is now an oriented manifold with an internal, 2-sphere boundary, whose topological 
complications are again confined to a compact region. Space-time is bounded to the future 
and past by two (partial Cauchy) surfaces M ± , extending to spatial infinity (see Figure 1). 

Following Definition 2 of weakly isolated horizons, we will equip A with a fixed equivalence 
class of vector fields [£] which are transversal to its 2-sphere cross-sections (where, as before, 
£ ~ £' if and only if £ = c£' for a constant c). It is also convenient to fix an internal null 
tetrad {£} , n 1 , m 1 , m 1 ) on A, each element of which is annihilated by the fiducial, flat internal 
connection d. 
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FIG. 1. The region of space-time M under consideration has an internal boundary A and is 
bounded by two partial Cauchy surfaces M which intersect A in the 2-spheres S' ± and extend to 
spatial infinity i°. 

The permissible histories consist of smooth triplets (e, A, A) on Ai satisfying boundary 
conditions at infinity and on A. The boundary conditions at infinity are, as before, the 
standard ones which ensure asymptotic flatness. Since the asymptotic behavior and bound- 
ary integrals at infinity play only a secondary role in our analysis, we shall not spell out 
the precise fall-off requirements. At A, the histories are subject to three conditions: i) the 
tetrads e should be such that the vector field £ a := £ ! ej defined by each history belongs to 
the equivalence class [£] fixed on A; ii) the tetrad e and the gravitational connection A should 
be such that (A, [£]) is a weakly isolated horizon for the history; and, hi) the electromagnetic 
potential A is in a gauge adapted to the horizon, i.e., C^A = 0. 

Remark: In space-time, we have the freedom to perform a local, internal Lorentz rotation 
on the tetrad ej (and the gravitational connection A a i J ). All these tetrads define the same 
Lorentzian metric g a b- Since i a is required to be a null normal to A, the permissible gauge 
rotations are reduced on A to the sub-group (IR + x -E 2 )i oc of local null rotations preserving 
the null direction field £. (Here H + is the group of rescalings of £ a , n a which leaves m a fixed 
and E 2 is the 3-dimensional Euclidean group consisting of rotations in the £-m, £-m and 
m-m planes.) Condition i) above — dictated by the existence of a preferred equivalence 
class [£] in Definition 2 — further reduces the internal gauge freedom to IR + x (E 2 oc ), i.e., 
reduces the group IR^ C of local £-n rescalings to the group 1R + of global rescalings. Thus, 
while any one space-time {A4,g a b), still defines infinitely many histories due to the freedom 
of tetrad-rotations, this freedom is somewhat reduced at A because of the structure fixed by 
the boundary conditions.^ 



6 Nonetheless, from a space-time perspective, the multiplicity of histories can still be rather sur- 
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Given any tetrad e", the internal null vectors (JL 1 , n 1 , m 7 , m 7 ) fixed on A trivially provide 
a null tetrad (£ a ,n a ,m a ,rn a ). In terms of these vectors, we can express E /J as: 



2£ [ V ] 2 e + 2n A (im£ [/ m J] - im£ [/ m J] ), (III. 2) 



where, as before, \ = im Am is the pull-back to A of the natural alternating tensor on 
the 2-sphere S of integral curves of £ a associated with the given history. The weak isolation 
of (A, [£]) restricts the form of the connection A at A. To see this, recall that one of the 
equations of motion requires the connection D defined by A to have the same action on 
internal indices as V. Hence, V a £i = d a £i + A a i J £j = A a i J £j where, in the second step 
we have used the fact that the flat derivative operator, d, has been chosen to annihilate the 
internal tetrad on A. Since V a e£ = by definition of V, and V ' a £ b = ui a £ b (see [11.4[ ) it follows 

that Ajj£ j = u£j. Hence, on A, A has the form: 

A u = -2l { jnj\u + C u , (III.3) 

where the 1-form Cu satisfies Cij£ j = 0. 

With this background material at hand, we are now ready to consider variations of the 
action ( III.1|) in the presence of an inner boundary representing a weakly isolated horizon. 



A key question is whether a new surface term at the horizon is necessary to make the 
variational principle well-defined. We will show that, thanks to the zeroth law, such a term 
is not needed. 

A simple calculation yields: 

8S(e, A, A) = / Equations of Motion • 5<p l — [ S /J A 5A L j - — [ 5 A A *F. (III.4) 

JM 16llG J A 47T J A 

where, in the first term, stands for the basic fields (e, A, A) in the action. Note that, as 
in the case without an internal boundary, the variation of the boundary term at infinity 
precisely cancels the contribution arising from the variation of the bulk terms. 

In order to show that the action principle is viable, it is necessary to show that the terms 
at the horizon vanish due to the boundary conditions imposed there. Let us begin with the 
gravitational term. Using ( |111.2j ) and ( |111.3| ) it can be re-expressed as 



1 



InG J a 



5co A \. (III.5) 



Since \ is the pull-back to A of the alternating tensor on the 2-sphere S of integral curves of 
£, it follows that Ce 2 e = 0. Furthermore, the weak isolation of the horizon ensures Cgu = 



prising. For example, if g a \, is the Schwarzschild metric with mass M > 0, there is a history in 
which the surface gravity Km is positive and another in which it is zero. This redundancy can be 
eliminated by working with isolated, rather than weakly isolated horizons. 
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and, since the null normal £ a defined by any tetrad belongs to the fixed equivalence class [£} 
at the horizon, we have 5£ = c$£ for some constant c$. These two facts imply Cgbuj = 0. Thus 
the entire integrand is Lie dragged by £. In the variational principle, however, all fields are 
fixed on the initial and final hypersurfaces, say M^. In particular, Sou necessarily vanishes on 
the initial and final cross sections of the horizon. Therefore, the integrand in ( |111.5| ) vanishes 
on the initial and final cross sections and is Lie dragged by £. This immediately implies 
( |111.5| ) is zero. 

Let us now consider the electromagnetic term. Since every A is in a gauge adapted to the 
isolated horizon, C e A = 0. Furthermore, 5£ a = cs£ a , so we conclude £i5A = 0. Next, ( |II.20| ) 
ensures Ce*F = 0. Thus, the integrand of the electromagnetic surface term is Lie dragged 
by £ a . An identical argument to the one presented above implies that the electromagnetic 
surface term in (|111.4| ) also vanishes. Therefore, the variation of the action ( |111.1| ) continues to 
yield Einstein-Maxwell equations in spite of the presence of an inner boundary representing 
a weakly isolated horizon. 

It is instructive to re-examine the key step in the above argument. Suppose we only had a 
non-expanding horizon. Then, the gravitational surface term could still be reduced to ( |111.5| ), 
and 2 e and *F would still be Lie-dragged by £. However, in this case, we could not argue 
that uj and A are also Lie-dragged. As we saw in Sections |11 L| and |11 C , these conditions are 
equivalent, respectively, to the constancy of the surface gravity Km and the electromagnetic 
potential on A. In this sense, given a non-expanding horizon as the inner boundary, 
the gravitational and electromagnetic zeroth laws are the necessary and sufficient conditions 
one must impose for the viability of the standard, first order, tetrad action principle. 

Remark: Note that is not the unique viable action for the problem: as usual, there 

is freedom to add suitable boundary terms without affecting the viability. Specifically, we 
are free to add any horizon boundary term which is composed entirely of fields which are 
Lie dragged by £, for example the intrinsic horizon metric q a b and fields u, 2 e, and A. Then, 
due to the argument given above, the new action would also be viable. However, as usual, 
this freedom will not affect the definition of the symplectic structure which underlies the 
Hamiltonian treatment of the next section. 



IV. CO VARIANT PHASE SPACE 

Let us now construct the phase space of space-times containing weakly isolated hori- 
zons. In the next section, we will use this framework to construct Hamiltonians generating 
suitable time translations and define the energy of an isolated horizon. In ||, the phase 
space was constructed by performing a Legendre transform of the action. This procedure 
leads to a 'canonical' framework in which the phase space consists of configuration and mo- 
mentum variables defined on a spatial hypersurface. With the self-dual connections used in 
||, the gravitational configuration variable turns out to be a connection and its conjugate 
momentum, a 2-form so that the Hamiltonian description can again be given in terms of 
forms. With the full Lorentz connections now under consideration, the situation turns out 
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to be more complicated. Specifically, one encounters certain second class constraints and, 
when these are solved, one ends up with the same canonical phase space that one would 
have obtained through a second order formalism. In the Hamiltonian framework, then, the 
simplicity we encountered in Section |T| is lost. More specifically, constraint functions and 
Hamiltonians now contain terms involving second derivatives of the basic canonical vari- 



ables which make variations rather complicated. (For details, see chapters 3 and 4 in |21 



Therefore, in this section we will not use a Legendre transform. Instead, we will construct 
the 'covariant phase space' from the space of solutions to field equations (see, e.g., 0]). As 
in Section |TJ, all expressions will now involve only the basic form-fields and their exterior 
derivatives and variational calculations will continue to be simple. 

To specify the phase space, let us begin as in Section |T| by fixing a manifold M. with 
an internal boundary A (see figure 1). As before, we will equip A with an equivalence 
class [£] of vector fields transverse to its 2-sphere cross-sections. To evaluate the symplectic 
structure and Hamiltonians, we will often use a partial Cauchy surface M in the interior of 
M. which intersects A in a 2-sphere S. Points of the covariant phase space V will consist of 
histories considered in Section |TJ which satisfy field equations. More explicitly, V consists 
of asymptotically flat solutions (e, A, A) to the field equations on M. such that i) the vector 
field £ a := £ 7 eJ belongs to the equivalence class [£] fixed on A, ii) in each solution, (A, [£]) is 
a weakly isolated horizon; and, iii) the electromagnetic potential A is in a gauge adapted to 
the horizon, i.e., CiA = 0. 

Our next task is to use the action ( |111.1| ) to define the symplectic structure Q on T. It is 
convenient to make a brief detour and first introduce two new fields which can be regarded 
as 'potentials' for the surface gravity Km and the electric potential Given any point 

(e, A, A) in the phase space T, let us define scalar fields ip and x on A as follows: 

i) Cnjj = (£ ■ uj) = Km and Cpx — {£ • A) = —Qm] and 

ii) ip and x vanish on S~~ , the intersection of M~ with A|] 

Note that these conditions associate with each point of V a unique pair (if), x) on A and in 
the 'extremal' case Km = 0, ip vanishes identically. 

We wish to use the standard procedure ]22] involving second variations of the action to 
define the symplectic structure.^ Let us recall the main steps of this procedure. One first 
constructs the symplectic current J: Given a point 7 in the phase space V and two tangent 



7 Condition ii) serves only to fix the freedom to add constants to ip an d X- One could envisage 
replacing it by a different condition. Our results will be insensitive to this choice. 

8 Actually this procedure provides a pre-symplectic structure, i.e., a closed 2- form on the phase- 
space which, however, is generally degenerate. The vectors in its kernel represent infinitesimal 
'gauge transformations'. The physical phase space is obtained by quotienting the space of solutions 
by gauge transformations and inherits a true symplectic structure from the pre-symplectic structure 
on the space of solutions. The 2-form f2 introduced below is indeed degenerate. However, for 
simplicity, we will abuse the notation somewhat and refer to £1 as the symplectic structure. 
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vectors 5% and 6 2 at that point, J provides a closed 3-form 7(7; 61, 6 2 ) on M.. Integrating dJ 
over the part M. of space-time under consideration, one obtains 



= /_ dJ( T ,5 h 8 2 ) = I _ J. 

JM JdM 



Now, if there is no internal boundary, one can choose M. to be a region bounded by any two 
Cauchy surfaces Mi and M 2 so that the boundary is given by dM. = Mi U M 2 U Too, where 
Too is the time-like 'cylinder at infinity'. In simple cases, the asymptotic conditions ensure 
that the integral J T J(7;5i,5 2 ) vanishes. Then, taking orientations into account, it follows 
that J M 7(7; 61,62) is independent of the choice of Cauchy surface M. One then sets the 
symplectic structure to be 



n\ 1 (s 1 ,6 2 ) = I J(r,6u6 2 ). 

JM 



In our case, the second variation of the action ( |III.1| ) yields the following symplectic 
current: 

J(T, 61, 62) = -^[5i£ /J A 6 2 Ajj - 6 2 Y, IJ A SiAu] - ^-[5i*F A 6 2 A - 5 2 *F A ^A]. (IV. 1) 

107rG 47T 

Using the fact that the fields 7 = (e, A, A) satisfy the field equations and the tangent vectors 
6\, 6 2 satisfy the linearized equations off 7, one can directly verify that J(-y; 61, 6 2 ) is in fact 
closed as guaranteed by the general procedure involving second variations. It is now natural 
to choose M. to be a part of our space-time M. bounded by two partial Cauchy surfaces 
Mi, M 2 , the time-like cylinder r ra and a part A of the isolated horizon bounded by Mi and 
M 2 . Again, the asymptotic conditions ensure that the integral of J over vanishes. Hence, 



/ +/ + L)J(r,Si,s 2 ) = 0. 

J Mi JM 2 JA 



'M 2 

However, this does not immediately provide us the conserved symplectic structure because 
the integral of J over A does not vanish in general. Since the isolation of the horizon implies 
that there are no fluxes of physical quantities across A, one might expect that, although 
non-zero, the integral over A would be 'controllable'. This is indeed the case. Using the 
forms ( |IH.2|) and (|III.3|) of S and A on A and the definitions of the potentials ip an d Xi ft 1S 



easy to verify that the pull-back of the symplectic current to A is itself exact: 

J(j;Si,S 2 ) = dj(y,5i,6 2 ) 
where the 2-form j on A is given by: 

J (7, *i, $0 = ^( 2 e) - 62^) } + ^[6 lX 6 2 *F - 62x61 *F] 

07TG 47T 

Hence, if Mi and M 2 intersect A in 2-spheres 5*i and S2 respectively, we have: 
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J(r,Si,s 2 ) = -( + )j(r,Sx,5 2 ). 

J5i JS 2 

The negative sign appearing in the above expression is due to the choice of orientation of 
Sa, which is induced from M rather than from A. Using these results we can define the 
symplectic structure as: 

/ [6 1 E IJ A 5 2 A U - 5 2 Z IJ A S 1 A IJ ] + ^<f 5 2 ^ - 5 2 {\) 5^} 

lb7rG JM 07TG Js 

-— I [6i*F A 5 2 A - 5 2 *F A 5iA] + — / [6SF 5 2X - 5 2 *F 8 lX ] (IV.2) 

47T JM 47T JS 

Again, using field equations one can directly verify that the right side of (|IV.2j ) is independent 
of the choice of the partial Cauchy surface; the symplectic structure is 'conserved'. We will 
use (r, fl) as our covariant phase space. 

Note that, even though the action did not contain a surface term at the horizon, the sym- 
plectic structure does. So, the overall situation is the same as in the undistorted, non-rotating 
case considered in M. Finally, our discussion of the action principle and our construction 
of the covariant phase space is applicable to all weakly isolated horizons A; nowhere did we 
have to restrict ourselves to the non-rotating case. 



V. HAMILTON! AN EVOLUTION AND THE FIRST LAW 

To discuss the first law, we must first define horizon energy, which in turn requires a time 
evolution field t a on j\A. Given a vector field t a satisfying appropriate boundary conditions, 
St := (C t e, C-tA, £tA) satisfies the linearized equations for any 7 := (e, A, A) in Y and thus 
defines a vector field on This 8 t can be interpreted as the infinitesimal generator of time 
evolution on the covariant phase space. It is then natural to ask if this vector field is a 
phase space symmetry, i.e., if £ s t Q vanishes everywhere on T. The necessary and sufficient 
condition for this to happen is that there exist a function H t — the Hamiltonian generating 
the t-evolution — such that 

6H t = n(6,5 t ) (V.l) 

for all vector fields S to V. On general grounds, one expects H t to contain a surface term 
-^adm a t infinity representing the total (i.e. ADM) energy, and a surface term E A at the 
horizon which can be interpreted as the horizon energy, both tied to the evolution field t a . 



9 In the Lie-derivatives, the internal indices are treated as scalars; thus Ct^i = t b db e a + e b^a^ fc - To 
make St a well-defined vector field on T, we now exclude from j\4. and let Ai and A be without 
future and past boundaries. Whenever needed, these boundaries, M ± and S^, can be added by 
taking the obvious closure of J\4 . 
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A key question then is to specify the appropriate boundary conditions on t a . It is clear 
that, at infinity, t a should be an asymptotic time-translation, i.e., should approach a time- 
translation Killing field of the flat metric used to specify the boundary conditions. At the 
horizon, on the other hand, the metric is not universal and the space-time defined by a generic 
point 7 of the covariant phase space does not admit any Killing field near A. Therefore, 
specification of the boundary conditions at A is not as straightforward as that at infinity. It 
is for this reason that we now assume that (A, [£}) is a non-rotating, weakly isolated horizon 
for all points 7 = (e, A, A) of the phase space. The problem of specifying the appropriate 
boundary conditions on t a in the rotating case is more complicated . However, it has been 
addressed successfully and will be discussed in [|13[ . 



Recall that the internal boundary A of M. is equipped with a specific equivalence class 
[£] of vector fields. As discussed in Section [TI|, these £ are the isolated horizon analogs of 
constant multiples of Killing fields on the Killing horizons in static space-times. Therefore, 
in the non-rotating case, it is natural to demand that, on A, the evolution vector field t a 
should belong to the equivalence class [£] . This automatically ensures that (£ t e, C t A, C t A) 
satisfy the appropriate boundary conditions to define a tangent vector at each point of the 
phase space V . However, unlike at infinity, the geometry at the horizon is not fixed once 
and for all. Therefore, it is natural to allow the precise value of the evolution vector field 
t a on A to vary from one point of the phase space to another. In more familiar terms, 
this corresponds to allowing the (boundary values of) lapse and shift fields to depend on 
dynamical fields (e, A, A) themselves, a procedure routinely used in numerical relativity and 
gauge-fixed calculations in canonical gravity. Following the current terminology in numerical 
relativity, we will refer to such t a as live evolution vector fields. The use of live fields turns 
out to be necessary to ensure that 5 t is a phase space symmetry, i.e., yields a Hamiltonian 
evolution on (T,J1). 

Let us fix a live t a whose restriction to the horizon belongs to the equivalence class [£} at 
all points of the phase space. To analyze if 5 t is a Hamiltonian vector field, it is simplest to 
compute the 1-form X t on V defined by 

X t (5) = Q(5,6 t ). (V.2) 

Now 5t is Hamiltonian — i.e., £ $ t Q = on T — if and only if X t is closed, i.e., 

<flX t = 

where «1 denotes the exterior derivative on (the infinite dimensional) phase space T. If this 
is the case then, up to an additive constant, the Hamiltonian is given by 

dH t = X t . 



To calculate the right side of ( |V.2|) , it is useful to note the following identities from 
differential geometry: 

£ t A = t ■ F + D(t ■ A) C t E = t- DE + D(t-E) -\(t-A),E] 

(V.3) 

C t A = t -F + d(t ■ A) C t *F = t • (d*F) +d(*-*F) 
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Using these, the field equations satisfied by (e, A, A) and the linearized field equations for 
5, we obtain the required expression of X t \ 



X t (5):=Q(5,Ct) (V.4) 
= — l -f Tr[(t.A)5Z-(t.Z)A5A}-±- f (f.A)*(*F)-(f.*F)A*A. (V.5) 

l07TG JdM 47T JdM 

Note that the expression involves integrals only over the 2-sphere boundaries Scx> and Sa of 
M, the partial Cauchy surface used in the evaluation of the symplectic structure; there is no 
volume term. 

The integrals at infinity can be evaluated easily by making use of the fall-off conditions. 
As one would expect, the electromagnetic term vanishes (because A falls off at least as 1/r 
while F falls off as 1/r 2 ) while the gravitational term yields precisely the ADM energy -Eadm 
associated with the asymptotic time-translation defined by t a . At the horizon, we can use 
equations (|III.2|) and ( 111.31) to show that t-T, contracted on internal indices with 8 A vanishes 



and ( |11.17| ) implies t ■ *F = leaving 

M5) = -7^[ (t ■ u)8( 2 e) - i / (t-A)8(*F)+8E t ABM 
87rG Jsa 47T Js A 

= -^G K W 5a * ~ $ W^a + 8E l ADM (V.6) 



where, in the last step, we have used the fact that both t ■ uo = km and t ■ A = — are 
constant on the horizon and the definition ( [II.19D of electric charge. 



The necessary and sufficient condition for the existence of a Hamiltonian is that X t be 
closed. Clearly, this is equivalent to 

— — dk (t) A <ffla A + d$ (t) A dQ A = , (V.7) 

where/Y\denotes the antisymmetric tensor product on T. Now, ( |V.7| ) trivially implies that the 
surface gravity Km and the electric potential at the horizon defined by t a can depend 
only upon the area and charge of the horizon. Other factors, such as the 'shape' of the 
distorted horizon, can not affect the values of km or Finally, ( |V.7|) is the necessary and 



sufficient condition that there exists a function E A , also only of a A and Qa such that 

SEi = ^ K {t) 8a A + $ (t) 8Q A . (V.8) 

Since E A is a function only of a A and Q A , it is a function of fields defined locally at the 
horizon. As noted before, it is natural to interpret E A as the horizon energy defined by the 
time translation t a . The total Hamiltonian is given by: 

H t = Ei DM - El (V.9) 

Let us summarize. Eq ( |V.8| ) is a necessary and sufficient condition for the 1-form X t to 
be closed. Therefore, the vector field 8 t on T defined by the space-time evolution field t a is 
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Hamiltonian if and only if the first law ( |V.8| ) holds. Thus ( |V.8| ) is a restriction on the choice 
of the live vector field t a : While any t a (which preserves the boundary conditions) defines 
an evolution flow on the phase space, it is only when 

1 

— — K W da A + $ (t )dQ A 

o7T(jr 

is an exact 1-form on V that this flow is Hamiltonian (i.e., preserves the symplectic struc- 
ture). At first, this restriction seems somewhat surprising because, in absence of internal 
boundaries, every vector field t a (which tends to a fixed Killing field of the flat metric at 
infinity) defines a Hamiltonian evolution. However, even in this context, there is no a priori 
reason to expect this tight correspondence to hold if one allows general, live vector fields 
t a whose boundary values at infinity can change from one space-time to another. Finally, 
we will see in Section |V11| that every space-time belonging to the phase space T admits an 
infinite family of vector fields t a for which X t is closed. Therefore, in particular, the first 
law does not restrict the 'background' space-times (or the variations 5) in any way. Indeed, 
for any space-time in our phase space, there is an infinite family of first laws, one associated 
with each permissible t a . 

We will conclude this section with a few remarks. 

i) Form of H t : The Hamiltonian ( |V.9| ) contains only surface terms. This may seem 
surprising because, in the canonical framework, the familiar Hamiltonian contains a volume 
integral consisting of a linear combination of constraints. While the volume term vanishes 'on 
shell' and does not contribute to numerical value of the canonical Hamiltonian on physical 
states, it is nonetheless crucial for obtaining the correct evolution equations since derivatives 
of the Hamiltonian transverse to the constraint surface are needed to construct the Hamil- 
tonian vector field. The covariant phase space, by contrast, consists only of solutions to the 
field equations whence the issue of taking 'off shell' derivatives never arises. In diffeomor- 
phism invariant theories, the Hamiltonian on the covariant phase space is always made of 
surface terms.Q If space-time has several asymptotic regions, the boundary term in each 
region defines the standard energy corresponding to that region. Therefore, in the present 
case, it was natural to interpret i? A as the horizon-energy defined by the t a evolution. Fi- 
nally, we should emphasize that we used a covariant phase space only for simplicity. The 
final results go through (and, in fact, were first obtained) in a canonical framework as well. 

ii) Comparisons: As noted in the Introduction, all treatments of the first law for non- 
rotating but possibly distorted horizons available in the literature refer to static space-times. 
The isolated horizon framework, by contrast, does not refer to a Killing field at all and thus 
allows a significantly larger class of physically interesting situations. On the other hand, 



In particular, therefore, Hamiltonians generating diffeomorphism which have support away from 
the boundaries vanish identically. Unlike their counterparts on the canonical phase space, the 
infinitesimal phase-space motions induced by such space-time vector fields are in the kernel of the 
covariant symplectic structure. 
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since it relies on a Hamiltonian framework, we cannot incorporate phenomenological matter 
if it does not admit a phase space description. Other treatments based on Hamiltonian 
methods generally restrict themselves to static space-times with a non-zero surface gravity. 
This assumption is essential there because those treatments use 'bifurcate' surfaces in an 
important way and these do not exist in the extremal static solutions where the surface 
gravity vanishes. In contrast, the results of this section, do not refer to a bifurcate surface 
and go through irrespective of whether km is non-zero or zero. In a realistic collapse, the 
physical space-time is not expected to have the bifurcate surface. The present analysis uses 
only the portion of the physical space-time in which the horizon has settled down with 
no further in-going radiation, rather than an analytical continuation of the near horizon 
geometry used in certain approaches. Finally, in contrast to other treatments, we have an 
infinite family of first laws, one for each evolution field t a for which St is a phase space 
symmetry. 

iii) Non-uniqueness of energy: Each permissible, live t a defines a horizon energy E\. At 
first it seems surprising that there is so much freedom in the notion of energy. Let us compare 
the situation at X + , which, like A, is null. There, we only have a 3-parameter freedom which, 
furthermore, can be eliminated simply by fixing a rest frame. How does this difference arise? 
Recall that energy is (the numerical value of) the generator of an unit time translation. 
At infinity, all 4-metrics in the phase space approach the same flat metric. Hence, we can 
simply fix a unit time-translation Killing field £" of that flat metric near infinity and use its 
restriction to X + as the unit time-translation for all metrics in the phase space. By contrast, 
there is no fixed 4-metric near A to which all the metrics in our phase space approach. 
Hence, we do not have the analog of only the equivalence class [£] is now common to 
all the metrics. If, for a given metric g a b in our collection, we select the time translation 
represented by a specific l a in [£], a priori we do not know which vector field l a in [£] would 
represent the 'same' time-translation for another geometry g a f,. One might imagine using the 
seemingly simplest strategy: just fix a l a a in [£] and demand that t a approach that £q for all 
points 7 in the phase space. Unfortunately, the strategy is not viable because such a t a fails 
to define a Hamiltonian evolution in Finally, if we restrict ourselves to globally static 
space-times, we can overcome this difficulty by always working with the Killing field which 
is unit at infinity. However, in absence of global Killing fields, the behavior of the evolution 
vector field t a near the horizon is unrelated to its behavior near infinity. Nonetheless, as we 
shall show in Section [VI]] , if one has sufficient control on the space of static solutions of the 
theory under consideration, it is possible to select a preferred energy function on the phase 
space and use it as the mass of the isolated horizon. In the Einstein-Maxwell case, all static 
solutions with horizons are explicitly known whence the strategy is viable. 



n As we saw above, a necessary condition for St to be a Hamiltonian vector field on T is that km 
is a function only of a a and Therefore, if we can find a tangent vector S in the phase space 
with SaA = SQa = but Sn^ ^ 0, t a can not define a Hamiltonian evolution. It is easy to find 
such a tangent vector 5 for this t a . 
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VI. YANG-MILLS FIELD 



In the previous three sections we restricted our attention to Einstein-Maxwell theory. 
We will now indicate how Yang-Mills fields can be included. This section is divided into 
three parts. In the first, we discuss restrictions on the Yang-Mills fields due to the horizon 
boundary conditions and introduce the notion of a Yang-Mills gauge adapted to the horizon'. 
In the second part, we discuss the action principle and construct the covariant phase space 
for Einstein- Yang-Mills theory. Using this formalism, in the third subsection, we introduce 
a Hamiltonian generating time evolution and extend the first law to the Yang-Mills case. 



A. Preliminaries 

We will restrict ourselves to compact gauge groups G and Yang-Mills connections defined 
on trivial bundles. Since the bundle is trivial, the connection gives rise to a smooth, globally 
defined Lie algebra valued 1-form, A. As usual, the Yang-Mills derivative operator D will be 
defined as DA = d\ + [A, A], where d is a flat Yang-Mills connection, and the field strength, 
F, via 

F := dA + A A A. (VI. 1) 

The stress energy tensor, T, is given in terms of the field strength as 

T ab = -^-[F ac l F b c t - ^FjF^l (VI.2) 

47T 4 

where the label i runs over the internal indices in the Lie algebra of the group G. 

Let us begin by examining how the isolated horizon boundary conditions restrict the 
form of the field strength, F, on A. Since the Yang-Mills stress energy tensor has the same 
form as the Maxwell one, ( |IL14| ), the analysis is completely analogous to that of Section 



11 C| . Therefore, we shall not include derivations of the results, but instead highlight the 
differences. 

Recall that, on a non-expanding horizon, R ab £ a £ b = whence T ab £ a £ b = 0. This has 
several consequences for the Yang-Mills field. In particular, one concludes 

£ a F ab = and £ a *F ab = 0. (VI.3) 

These two restrictions guarantee there is no flux of Yang-Mills field across the horizon. 
Making use of the specific form of the stress energy tensor, we also conclude 

$02 = ^R ab m a m b = and $ 20 = ^R ab rn a m b = 0. (VI.4) 

Our next task is to define the Yang-Mills equivalents of the electric and magnetic charges 
of the horizon. Naively, one might consider integrating F and *F over a 2-sphere cross section 
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of the horizon as was done in the Maxwell theory. However, these 2-forms now have a free 
internal index and are only gauge covariant rather than gauge invariant. Since there is no 
preferred internal basis at the horizon, the integrals would fail to be well-defined. Therefore, 
we must look for 2-forms which are gauge invariant. A natural quantity to consider is 
the norm of F, defined by the Killing-Cartan form Kij on the Lie-algebra of G and the 
(contravariant) \ on the horizon: 



\ ■ F)*( 2 e • F) j Kij\ 2 (VI.5) 

(Although the contravariant \ is ambiguous up to terms of the type £^- a V b ^ where V a is any 
vector field tangential to A, this ambiguity does not affect |F| because of ([VI. 3| ). The norm 
of *F is defined analogously. These two quantities are gauge invariant and allow us to define 
the electric and magnetic Yang-Mills charges of the horizon: 

QYM.^_J_I |* F | 2 e pYM.^l I | F |2 £ (VL6) 

A 4vr Is a 4tt Js a 

Recall that the unusual signs in the definitions of the charges arise due to the orientation of 
the S/\ — the normal to the two sphere is inward pointing. In Maxwell theory, the magnetic 
charge is zero unless we consider either connections on non-trivial bundles or allow 'wire 
singularities'. As is well known, this is not true for Yang-Mills theory: the magnetic charge 
can be non-zero even if we restrict attention to smooth fields on a trivial bundle. 



We would now like to verify that the charges defined in ( [VI. 6|) are independent of the 
cross section of the horizon on which the integration is performed. The isolated horizon 
boundary conditions guarantee this is the case. First, recall the geometric identity 

C e F = £ ■ DF - [(£ ■ A), F1 + D(l ■ F) . (VI.7) 

A similar expression for *F is also true. The first term on the right hand side vanishes due 
to the field equations and the third term is zero due to the previous restriction on F, ( VI. 3|) . 
Therefore at the isolated horizon, 

£ e F = -[{£■ A), F] and C /F = -[(£■ A) ,*F] . (VI.8) 

In the Maxwell case, F and *F are Lie dragged by i. However, for non-Abelian fields, this is 
not a gauge invariant statement; the terms on the right hand sides of ( |VI.8| ) are necessary 
for gauge invariance. Although the field strength and its dual are not Lie dragged along £ a , 
recalling that Ce 2 e = and using the cyclic property of the trace, it is straightforward to 
demonstrate that their norms are: 



£*|F|=0 and £ e \*F\ = 0. (VI.9) 

This result, along with ( |V1.3| ), guarantees that the charges Q\ M and P™ are independent 
of the choice of cross section Sa of the horizon. 
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Let us consider the remaining components of the Yang-Mills field. The boundary con- 
ditions place no restrictions on F a i,n a m b and *F a i,n a m b at all. As in the electromagnetic 
case, these components describe the radiation flowing along the horizon. The isolated hori- 
zon boundary conditions allow radiation arbitrarily close to — and even at — the horizon, 
provided none crosses it. 

We have so far restricted our attention to the field strength and its dual. However, in 
the action principle and phase space, the basic variable will be the Yang-Mills connection 
A. Let us begin with the definition of the Yang-Mills equivalent of the electric potential. 
Recall that, given an £, the electric potential was defined in Section |II Q as $m = — (£ • A). 
This definition is not appropriate in the Yang-Mills case since the resulting potential has a 
free internal index and is therefore not gauge invariant. Instead, we define the Yang- Mills 
potential, $™ to be negative the norm of (£ ■ A) : 

$™:=-\(£-A)\. (VI.10) 

This gives us a gauge invariant potential at the horizon. 

As in Maxwell theory, we need to constrain the form of A at the horizon. Several 
considerations motivate our choice of these boundary conditions. First they must be chosen 
so that the action principle is well defined. Second, if the gauge group is U(l) the boundary 
conditions should reduce to those given in Section [11 G for the electromagnetic field. Finally, 
we should be able to show that the Yang-Mills electric potential is constant on the horizon. 
These considerations suggest the following definition: 

Definition 6: The connection A will be said to be in a gauge adapted to the isolated horizon 
(A, [£]) if it satisfies the following two conditions 

• (i) The Yang-Mills potential is constant on the horizon 

t^f = (VI. 11) 



(ii) The dual of the field strength (*F) and [I ■ A) point in the same Lie algebra 
direction, 

{£ ■ A) { oc ( 2 e • *F)* (VI.12) 



on the horizon. 

These boundary conditions satisfy the requirements discussed above. First, it is straight- 
forward to show that in the U(l) case, condition (i) is equivalent to requiring £ eA = and 
(ii) is redundant. Second, as we shall see in Section |VI B| , these boundary conditions are also 
sufficient to make the variational principle well defined. 

It is not difficult to show that these conditions can always be satisfied. The remaining 
gauge freedom is simply A — > g^Ag + g^dg, where g satisfies C ig = 0. 
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B. Action and Phase Space 



In this section we will consider the first order action for Einstein- Yang-Mills theory 



on the manifold jW described in Section [III B| . The basic fields will consist of the triplet 
(e^, A a i J , A*), where e f a and A a j J are the tetrad and Lorentz connection and A^ is the Yang- 
Mills connection. The gravitational fields, and A a j J satisfy the same boundary conditions 



(at A and infinity) as in Section pIB|. Furthermore, we require the Yang-Mills fields to be 
in a gauge adapted to the horizon and assume they fall off sufficiently fast at infinity.^ The 
Einstein- Yang-Mills action is 

S{e, A, A) = f S /J A F u + — — f S /J A A u — / Tr[F A *F] . (VI.13) 



16nG Jm 16nG j 671 JM 

The gravitational part of the action has previously been discussed in Section |T|, therefore 
we shall only consider in detail the Yang-Mills terms and verify the variational principle is 
well defined. Taking into account the results of Section |HI| , a variation of the action can be 
expressed as 

6S(e, A, A) = f Equations of Motion • 8(j) - — [ Tr [5 A A *F]. (VI. 14) 

Jm 47r J a 

We must demonstrate the boundary term at A vanishes due to the conditions imposed on 
the Yang-Mills fields. Using ( |V1.3| ) and ( |V1.12| ), one can show that the trace in ( |V1.14| ) can 
be replaced by a product of norms: 

Tr[5A A *F] = - cs^f) | *F | A e , 

where A e = n A \ is the volume form on A and S£ = c$£. 

In the action principle, variations are performed keeping data fixed at the initial and final 
slices. In particular <5$™ and cs vanish there. However, the boundary conditions guarantee 
that and hence its variation, is constant on A. Since, 5<&™ and c$ vanish on the initial 
cross section of the horizon and are constant, it follows that = and cs = 0. Therefore, 

the Yang-Mills horizon boundary term vanishes; the action principle is well defined in the 
presence of Yang-Mills fields. As in the Einstein-Maxwell case, boundary conditions played 
a crucial role in demonstrating the viability of the action. 



12 More specifically we require the fall off conditions on the Yang-Mills connection to be such that 
all integrals, in particular the symplectic structure, are finite and yet the asymptotic electric and 
magnetic charges are not forced to vanish. While it is not trivial to meet these conditions (for 
examples the conditions used in [|5,12] appear not to lead to a well-defined symplectic structure) 
they can be met. However, as in the rest of the paper, for brevity, we will not spell out the boundary 
conditions at infinity in detail. 
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We now wish to construct the covariant phase space and symplectic structure. As be- 
fore, points in covariant phase space T will consist of histories which satisfy the Einstein- 
Yang-Mills field equations, appropriate falloff conditions at infinity and the isolated horizon 
boundary conditions at A. Before proceeding further, we shall once again need to introduce 
an additional field at the horizon. This can be regarded as 'potential' for the Yang-Mills 
potential $ YM . Given any point (e, A, A) in the phase space T, let us define the scalar field 
X on A as follows: i) Cex = — j an d ii) X vanishes on S~ , the intersection of M~ with 
A. These conditions are identical to those imposed in the Maxwell case. 

Once again, we take second variations of the action in order to obtain a symplectic 
structure. Since the gravitational terms are exactly the same in Section [TTT], we shall only 
describe in detail the Yang-Mills part of the symplectic structure. The second variation of 
the action ( [VI. 13|) yields the following symplectic current: 

J( T , <5i, S 2 ) = J grav - -^Trf^F A 5 2 A - 5 2 *F A MJ. (VI.15) 

Using the fact that the field equations and linearized field equations are satisfied, one can 
directly verify that J{^;5i,5 2 ) is indeed a closed 3-form. We again choose the spacetime 
region of interest, Ai, to be that part of the spacetime M. bounded by M\, M 2 , infinity and 
a portion A of the isolated horizon. Integrating dJ over J\A and using asymptotic falloff 
conditions, we obtain 



/ + / + L)J(T,Sqi,6 2 ) = 0. 



The integral of J over A does not vanish but, as in Section |TJ, the pull back of J to A is 
exact. Therefore, we can express the integral over A of J as integrals over the initial and 
final 2-spheres S\ and S 2 . Using these results, and keeping track of orientations, we obtain 
the symplectic structure: 

f2| 7 (5i, 5 2 ) = 

167rG Jm 87rG Js 

_J_ f [S^F A 5 2 A — 5 2 *F A 5iA] + ±- <f [8 t (|*F | \)5 2X - S 2 (|*F|^ lX ]. (VI.16) 

47T JM 47T JS 

Full use of the isolated horizon boundary conditions has been made in obtaining this sym- 
plectic structure. In particular, to obtain the given form of the Yang-Mills surface term, we 
have used the fact that (i ■ A) and ( 2 e • *F) point in the same direction in the Lie algebra. 
Using field equations one can directly verify that the right side of (|VT.16|) is independent of 



the choice of the partial Cauchy surface; the symplectic structure is 'conserved'. We will use 
(T, Q) as our covariant phase space. 
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C. Hamiltonian and First Law 



In this subsection we will generalize the arguments of Section to obtain an expression 
for the energy of the horizon in Einstein- Yang-Mills theory. To do so, we must specify a 
time evolution vector field t a . As before, we require t a to be a member of the preferred 
equivalence class [£ a ] at the horizon (this again requires restriction to non-rotating isolated 
horizons) and approach unit time translation asymptotically Given t a we can calculate the 
infinitesimal generator of time evolution, 5t = (£te, CtA,£ t A.), and determine whether it is 
Hamiltonian. Recall that 5t is Hamiltonian if and only if the 1-form X t on the phase space 
defined by 

X t (5) :=n(8,6 t ) (VI.17) 

is closed. Let us calculate X t . The gravitational part will be identical to the expression 
obtained in Section [V|, therefore we shall concentrate on the Yang-Mills terms. As with 
the Maxwell field, the Lie derivatives of the Yang-Mills fields can be re-expressed using the 
following identities: 

C t A = t-F + D(t • A) C t *F = t ■ (D*F) — [£■ A, *F] + ~D(t ■ *F) . (VI.18) 

Making use of these expressions, the field equations satisfied by (e, A, A) and the linearized 
field equations for 5, we obtain the required expression for X t : 

X t {5) = —^-f (t.A)5Z-(t.Z)A5A-±- f Tr[(t.A)CF)-(t.*F)A5A]. 

1D7TG JdM 47T JdM 

(VI.19) 

As before the expression involves integrals only over the 2-sphere boundaries Sea and 5*a 
of M; there is no volume term. The gravitational terms yield 5E t ADM at infinity and 
— (1/8ttG)k^)Sci^ at the horizon. The Yang Mills term at infinity vanishes due to fall- 
off conditions, therefore we need only calculate the Yang-Mills contribution at the horizon. 



This is composed of two terms, the second of which vanishes due to the restriction (|VI.3| ) 
which guarantees £ ■ ^F = 0. Since we are in a gauge adapted to the horizon, (£ ■ A) and *F 
point in the same internal direction. This allows us to replace the trace in the first term by 
norms: 

Tr[(t-A)5rF)] = -$^(rF| 2 e) 

and guarantees that is constant on the horizon. Making use of the definition of Yang- 
Mills electric charge, ( |VT.6| ), we obtain 

Xt{6) = ~8^G K{t)6aA - ^r^r + ^adm- 

Recall that the necessary and sufficient condition for the existence of a Hamiltonian is that 
X t be closed. Clearly, this is equivalent to 
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dK(t) A da A + <fl$^ A dQX M = 0, (VI.20) 



Once again, we conclude that the surface gravity K(t) and the Yang-Mills potential at 
the horizon defined by t a can depend cm/?/ upon the area a& and charge QX M °f the horizon. 
Finally, (|VI.20| ) is also the necessary and sufficient condition that there exist a function E l A , 
also only of cia and QX M sucn that 

^1 = K(t)*OA + ^Ql M . (VI.21) 

As before, i? A is interpreted as the horizon energy defined by the time translation t a . We 



conclude that the vector field 5t is Hamiltonian if and only if the first law, (|VI.21|) , holds. 
We will conclude this section with a few remarks. 

i) The derivation of the first law and its final form are completely analogous to those in 
the Einstein-Maxwell theory. By contrast, in the discussion of the first law for undistorted 
isolated horizons of ||12|| , certain restrictions were imposed on the permissible variations 5 in 
the Einstein- Yang-Mills case. In our treatment, subtleties arise only in the definition of a 
canonical horizon mass (see Section |V11 C| ) rather than the discussion of the first law itself. 

ii) Although the Yang-Mills magnetic charge P][ M will generically not be zero, no term 
involving 5P][ M arises in the first law. 



iii) How does our result compare with those previously available |23],|24|]? In |24| the first 
law for Yang-Mills fields is proved for globally stationary spacetimes and small perturbations 
from one such space-time to another. Assuming the Yang-Mills fields fall off sufficiently fast 
at infinity, (in the non-rotating case) the first law of |24| then reads 



SM = -^—K5a + / Tr[0 5*F] . 
8nG Juot 



Here, M is the ADM mass evaluated at infinity, while all terms on the right hand side are 
evaluated at the horizon. Due to a different gauge choice at the horizon, the authors define 
a Lie algebra valued potential and leave the term inside an integral. However, the 

general form of this first law is the same as ours. In this sense, our framework generalizes 



the results of 241 to non-static contexts. 



In ||23|| , the first law is proved for globally stationary spacetimes and arbitrary small 



departures therefrom. However, there are a number of important differences between these 



results and the ones obtained in this paper. In the non-rotating case, the first law of p3 
reads 

SM ADM + V5Q = ——jKda 

07TU 

where V and Q are the Yang-Mills potential and charge evaluated at infinity while k and a 
are of course evaluated at the horizon. Because of the non-Abelian nature of the Yang-Mills 
field, unlike in the Maxwell case, the charge Q evaluated at infinity is now quite different 
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from the charge evaluated at the horizon and, as in the Maxwell theory, the potential V 
evaluated at infinity has no direct bearing on the potential at the horizon. Furthermore, 
that calculation makes an essential use of the bifurcation 2-sphere and all fields are required 
to be smooth there. This restriction implies that the Yang-Mills potential at the horizon 
vanishes. (The same is true if one restricts the analysis of |23| to the Maxwell case.) 

The first law derived in the isolated horizon framework is valid also in presence of radi- 
ation in the exterior space-time region and makes no reference to the bifurcation 2-sphere. 
(Although we restricted ourselves to the non-rotating case, rotation has been incorporated in 
this framework in ||13|1 .) Furthermore, it has the aesthetically pleasing feature that all quan- 
tities that appear in ( |V1.21|) — including the energy E l A , the potential $™ and the charge 
Qa M — are evaluated at the horizon. In particular, one can now meaningfully consider 
the physical process version in which one does an experiment at the horizon by dropping a 
test particle/field and changing the horizon charge infinitesimally. More generally (|VI.21|) is 
genuinely a law governing the mechanics of the horizon. 



VII. HORIZON MASS 



For notational simplicity, we will say that a (live) vector field t a is permissible if it 
gives rise to a Hamiltonian evolution. We saw in Sections |V| and [VI| that each permissible 



vector field t a defines a horizon energy E f A . In the phase space framework, E A has a direct 
interpretation: it is the surface term at the horizon in the expression of the Hamiltonian 
generating the devolution. However, in many physical applications — such as the study of 
black hole mergers — one is interested in properties of a specific space-time, rather than the 
full phase space. Then, it is useful to have at one's disposal a canonical notion of energy, 
the analog of the ADM energy in the rest frame at infinity. This quantity could then be 
interpreted as the horizon mass. In this section, we will introduce this notion in detail. 
The discussion is divided into three parts. In the first, we consider the Einstein-Maxwell 
theory; in the second, we discuss dilatonic couplings 0; and, in the third, we analyze the 
Einstein- Yang-Mills system. 



A. Einstein-Maxwell theory 

In Section [V] we showed that t a is permissible if and only if ( | V . 7| ) holds on the phase 
space. We will now construct a large family of permissible evolution fields t a . Fix any 
regular function k of two variables gia and Qa- Then, given any point 7 = (e, A, A) of T, 
we define (the boundary value of) the vector field t a as follows. Consider the vector field i a 
on A defined by the tetrad, £ a = efi 1 , and denote by Km the surface gravity associated with 
it. Then, k q = cKm for some constant c. Let us set t a = c£ a . Repeating this procedure at 
each phase-space point 7, we obtain a live vector field t a with Km = k q . (The resulting c will 
be constant on A but a function on the phase space.) Next, consider the electro-magnetic 
potential, which is guaranteed to be constant on A by our boundary conditions but whose 
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value at any phase space point is so far completely free. We will now use ( |V.7| ) to fix it. 
Equation (|V.7|) implies: 

8Qa dctA 

Since = n Q is known, we can simply integrate the equation for $( t ) as a function of cza 
and Qa- Furthermore, the solution is unique if we impose the physical condition that 
should vanish whenever Qa = 0. Thus, starting from any regular function n Q of a a and Qa, 
we have obtained a permissible evolution field t a . Conversely, it is easy to verify that every 
permissible vector field arises via this construction. There is clearly a very large family of 
such live vector fields. 

An obvious question is if there is a 'canonical' or 'natural' choice of t a l We will now 
show that the answer is in the affirmative. Recall that, in the Einstein-Maxwell theory, 
there is precisely a 2-parameter family of globally static solutions admitting horizons: the 
Reissner-Nordstrom family. (Since A is required to be a globally defined connection on a 
trivial U(l) bundle, the magnetic charge is zero on the entire phase space.) Let us focus on 
this family. Denote by £ a the static Killing field which is unit at infinity. Its surface gravity 
is a specific function of a a and Qa- 

1 (, GQl\ 

k <« = 2^ y 1 - -rtJ ■ 

As before, Ra is the horizon radius, defined by a a = 47ri? A . We can therefore use in 
place of k in the above construction. The resulting permissible, live vector field agrees 
with £ a on the horizon of every static solution. This property is satisfied only if we set 

Next, we can 'integrate' ( |V.8| ) to obtain the horizon energy Efc. Although a priori there is 
the freedom to add a constant, we can fix it by requiring that the energy vanish as a a and Qa 
tend to zero. Indeed, we have no choice in this since one cannot construct a quantity with 
dimensions of mass from the fundamental constants that appear in the Einstein-Maxwell 
theory. (Einstein- Yang-Mills theory does admit such a constant and we will see in Section 
|VII C| that it leads to an interesting modification of the situation discussed here.) Let us 
define the horizon mass via 

M A = E\. 

To justify this definition, let us begin by restricting to static solutions. In each static solution, 
we are free to extend t a Q such that it coincides with the Killing field £ a . General symplectic 
arguments imply that, on any connected component of the space of static solutions, the 
numerical value of the total Hamiltonian, generating evolution along £ a , must be constant 
(see, e.g., [§].) In the Einstein-Maxwell case, there is a single connected component and, by 
the dimensional argument given above, the numerical value of the Hamiltonian must vanish 
on it. Hence, from ( |V.9|) it follows that, on any static solution, 
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H to = M ADM - M A = 0. 



On a general solution, of course, Madm would be greater than Ma, the difference being equal 
to the energy in radiation. If the horizon is complete in the future and time-like infinity i + 
satisfies certain regularity conditions, as in || one can argue that the difference is precisely 
the total energy radiated across X + and hence Ma equals the future limit of the Bondi mass. 
These considerations support our interpretation of Ma as the horizon mass. 

Finally, since we now have a canonical evolution field we can drop the suffix t on surface 
gravity and electromagnetic potential and write the first law (jV.8|) in the more familiar form: 

SM A = — — K5a A + $5Q A ■ 

07TG 

In contrast to treatments based on static space-times, the quantities that enter this law are 
all defined at the horizon. Therefore, as pointed out in , it is now possible to interpret this 
law also in the 'active' sense where one considers physical processes which increase the area 
and the charge of a given horizon. To our knowledge, the standard proofs of this physical 
version p,ID| are not applicable to processes in which the background has non-zero electric 



charge and the process changes it infinitesimally. 
We will conclude with a few remarks. 

i) In the above discussion, the permissible evolution field t a Q was constructed by setting 
t a Q = c£ a where c is given by k q = (1/2.Ra)(1 — G(Q A /R A ) 2 ) = ck^. For c to be well- 
defined, it is necessary that k vanishes whenever Km does. Therefore, for the mass to be 
well-defined, we must excise those points from the phase space at which Km vanishes but k 
does not. However, this is not a serious limitation. In particular, we still retain all static 
solutions including the extremal ones at which k q vanishes. 

ii) Since we have a specific k , we can use ( |V.7| ) to obtain the corresponding electro- 
static potential: $ Q = Q A /R A . Furthermore, by integrating ( |V.8|) it is easy to express Ma 
explicitly in terms of the horizon parameters: 

Mi ^ KOi + * ea ^( 1 + ^). (VII,) 

Thus, the functional dependence of Ma on the horizon parameters at any point of the phase 
space is the same as in static space-times. Note that this is a result of the framework, not 
an assumption. Its derivation involved two distinct steps. First, and most importantly, the 
first law ( |V.8| ) arose as a necessary and sufficient condition for the existence of a consistent 
Hamiltonian framework. Second, the freedom in t a was exploited in order to construct the 
preferred, permissible evolution field It is quite significant that Ma can be expressed so 
simply using just the parameters defined locally at the horizon even when there is radiation 
arbitrarily close to it. This fact is likely to play an important role in the problem of extracting 



physics in the strong field regimes from numerical simulations of black hole collisions JL8 
It is important to notice that although we made use of our knowledge of static solutions to 
arrive at a canonical t a Q and the mass function Ma, the final result ([VII.l ) makes no reference 
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to these solutions. is a simple function of the parameters which can be directly computed 
from the geometry of any one isolated horizon. 

iii) In the earlier work j8| on undistorted horizons, one restricted oneself to the preferred 
evolution field t a Q from the very beginning (although this vector field was selected using a 
different but equivalent procedure). Q The a priori freedom in the choice of a permissible t a 
was not discussed and the first law appeared only in the more familiar form, given above. 



B. Dilatonic coupling 

The Einstein-Maxwell-dilaton system was studied in some detail in the undistorted case 
in 0. We will revisit it here in the more general context considered in this paper because it 
brings out a subtlety in the definition of the horizon mass Ma and the associated first law. 

The dilaton is a scalar field <ft which can couple to the Maxwell field in a non-standard 
fashion. The coupling is governed by a constant a. If ct=0, one obtains the standard 
Einstein-Maxwell-Klein-Gordon theory and the situation then is completely analogous to 
the Einstein-Maxwell theory considered above. If a=l, the theory represents the low energy 
limit of string theory. In this case, there are some interesting differences from the Einstein- 
Maxwell theory considered in this paper. To bring out these differences, in this subsection 
we will set a—1. (The situation for a general value of a is discussed in || where one can 
also find details on the material summarized below.) 

In the standard formulation, the theory has three charges, all defined at infinity; the 
ADM mass Madm, the usual electric charge Qoo and another charge Qoo : 

Qoo = ^~ i *F and Q 00 = ^- I e~ 2 ^ *F . 

Q is conserved in space-time (i.e. its value does not change if the 2-sphere of integration is 
deformed) while Q is not. From the perspective of the isolated horizons, it is more useful to 
use a a, Qa, Qa as the basic charges^: 



13 This strategy seems to have generated a misunderstanding (see, e.g., [|29[j ) that the first law was 
obtained in |8|,p|,[T2| merely by identifying the parameters labeling a general isolated horizon with 
those of static horizons and then using the Smarr formulas available in the static context. This 
was not the case. Rather, static solutions were used only to select the appropriate normalization 
of the evolution vector field t a at the horizon. The Hamiltonian framework was then used to define 
the horizon mass without any reference to Smarr formulas. As in this section, the mass was then 
shown to reproduce the Smarr-type formulas on general horizons. 

14 In the undistorted case, the dilaton is constant on A and hence we can replace Qa by 4>a as in 
1- 
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Q A = — f *F, 

47T JS A 



and Qa 



47T JS A 



-20* 



F. 



Although the standard electric charge is not conserved in space-time, it is conserved along 
A whence Qa is well-defined. 

It is straightforward to extend the construction of the phase space to include the dilaton. 
The only difference is that the charge Q in equations ( |V.6| ) - ( |V.£| ) is replaced by Q. With 
this minor change, the discussion of the first part of Section |VII A| is also unaffected. Thus, 
given any function n Q of a a and Qa, we can construct a permissible, (live) evolution field t a . 

The difference arises in the next step where we constructed a preferred t a Q . With the 



dilatonic coupling, the theory has a unique j25[ three parameter family of static solutions 
which can be labeled by (a a, Qa, Qa)- As in the Reissner Nordstrom family, these solutions 
are spherically symmetric. In terms of these parameters, the surface gravity kiq of the static 
Killing field which is unit at infinity is given by: 
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QaQl 
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QaQa 



The problem in the construction of the preferred is that we need a function k which 
depends only on a a and Qa- Therefore, we can no longer set k = Krg) on the entire phase 
space because K(£\ depends on all three horizon parameters. 

To extract the mass function Ma on the phase space, we can proceed as follows. Let 
us foliate T by Qa = const surfaces. On each leaf, Kr^ trivially depends only on <2a and 
Qa and so we can set k = K^y Therefore, by the procedure outlined in Section |V11 A| , 
we obtain a (live) vector field t a Q and can define the mass Ma (7) = ^a(t) for all points 7 
on this leaf. Repeating this procedure for each leaf, we obtain a live vector field t a Q and a 
mass function Ma everywhere on T. However, the surface gravity K(t ) now depends on all 
three parameters, rather than just a a and Qa- Therefore, the first law (|V.8|) cannot hold 
for arbitrary variations 5 and consequently 5t„ fails to be a Hamiltonian vector field. Put 
differently, although there is a multitude of permissible, live vector fields, each leading to a 
first law, none of them can coincide with the Killing field £ a (which is unit at infinity) on all 
static solutions. This is a significant departure from the Einstein- Maxwell case considered 
above. 

Nonetheless, (modulo the caveat discussed in the first remark at the end of Section 
[VII A|) the above procedure does provide us with a well-defined mass function Ma on the 
entire phase space which can be expressed in terms of the horizon parameters as 



2G 



QaQl 



It equals Madm in static space-times and has other properties which motivated our inter- 
pretation of Ma as the horizon mass in the Einstein-Maxwell case. Since this function is 
well-defined on the entire phase space, we can simply vary it and express the result in terms 
of the horizon parameters. The result is: 
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SM A = — — k5cl a + $5Q A 

where k = Ku), $ 2 = (QaQa/R a ) an d Q\ = QaQa- Thus, although there is still a 
first law in terms of £" and M A , it does not have the canonical form ( |V.8| ) because is 



not a permissible vector field. More generally, in theories with multiple scalar fields ||26|| , if 
one focuses only on static sectors, one obtains similar 'non-standard' forms of the first law 
with work terms involving scalar fields. This reflects the fact that there is no permissible 
vector field t a , defined for all points of the phase space, which coincides with the properly 
normalized Killing field on all static solutions. In the undistorted case, the analysis was 
carried out only in terms of the vector field t a Q and the horizon mass M A @- The resulting 
first law had the above form. 

Alternatively, one can restrict oneself to variations 5 which are tangential to the leaves 
of the phase space foliation constructed above. Since t% is a permissible vector field for any 
one leaf, we obtain the standard first law 

5M A = (l/87rG)K ito) 5a A + $ (to) 5Q A 

for the restricted variations. The idea of using such restricted variations was suggested in 
|12|1 in the context of Yang-Mills fields (although the foliations and other details were not 
spelled out there). 

To summarize, because there is now a three parameter family of static solutions rather 
than two — or, more precisely, because the standard surface gravity in static space-times 
depends on a A , Q A and Q A — a canonical, permissible evolution field is no longer available. 
However, there is still a multitude of permissible evolution fields and corresponding first 
laws. Furthermore, one can still define a canonical mass function M A on the entire phase 
space. 



C. Yang-Mills fields 

In Einstein-Maxwell theory, with and without the dilaton, there is no way to construct 
a quantity with the dimensions of mass from the fundamental constants in the theory. The 
situation is different for Einstein- Yang-Mills theory because the coupling constant g has 
dimensions (LM) _1//2 . The existence of a quantity with units of mass has interesting conse- 
quences which we will now discuss. 

Let us begin with a summary of the known static solutions in Yang-Mills theory. First, 
the Reissner-Nordstrom family constitutes a continuous two parameter set of static solutions 
of the Einstein- Yang-Mills theory, labelled by (a A , QX M )- I* 1 addition, there is a 1-parameter 
family of 'embedded Abelian solutions' with (a fixed) magnetic charge P A , labelled by 
(a A ,P^). Finally, there are families of 'genuinely non- Abelian solutions'. For these, the 
analog of the Israel theorem for Einstein-Maxwell theory fails to hold |27j ; the theory admits 



static solutions which need not be spherically symmetric. In particular, an infinite family of 
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solutions labelled by two integers (%, n 2 ) is known to exist. All static, spherically symmetric 
solutions are known and they correspond to the infinite sub-family (ni,n 2 = 0), labelled by 
a single integer. However, the two parameter family is obtained using a specific ansatz, so 
there may well exist other static solutions. Although the available information on static 
solutions is quite rich, in contrast to the Einstein-Maxwell-Dilaton system, one is still rather 
far from having complete control of the static sector of the Einstein- Yang-Mills theory. 

The zeroth and first laws do hold in the Einstein- Yang- Mills case. At present, however, 
we can only hope to repeat the strategy used in the last two sub-sections to define a canonical 
mass function Ma on portions of the phase space. In order to define it on the full phase 



space, of the 'uniqueness' and 'completeness' conjectures of [12] will have to hold (possibly 
with a suitable modification.^ Nevertheless, new insight into the static solutions can be 
obtained by restricting attention to certain leaves of the phase space. The basic idea is 
taken from |12| but applied in a slightly different manner to the more general context of 



distorted horizons. 

Consider a connected component of the known static solutions, labelled by n = (711,712). 
This is a 1-dimensional sub-space of the phase space which we denote S^. Each point in Sft 
can be labelled by the value of the horizon area a&. Calculate the surface gravity for 
this family, where £ a is the static Killing field which it unit at infinity and set k = Ken in 
the construction sketched in Section [VII A| . We then obtain a live vector field t°, and the 
corresponding first law 

SE * = 8^G K ^ 6aA 

on the full phase space. 

When restricted to Sn, we can interpret E*£ as the horizon mass Ma and replace Ku\ 
by the function /^(oa) used in the literature: /3g = 2K( to )i?A- Then, by integrating the first 
law along S^, one obtains: 



2G jo 

where, we have used the fact that, since E^ is a surface integral at A, it vanishes as the 
horizon area goes to zero. Thus, the horizon mass is completely determined by f3ft(a&). 

Next, we use the fact that the Hamiltonian given by H to = E t ^ DM — E*£ (see ( |V.9| )) is 
constant on each connected, static sector if t„ coincides with the static Killing field on the 
entire sector. By construction, our t„ has this property. In the Einstein-Maxwell case, since 



5 For example it may be appropriate to restrict oneself to the class of space-times admitting 
isolated horizons which are complete in the future. Physically, this is the most interesting case since 
such horizons would represent the asymptotic geometry resulting from a gravitational collapse or 
black hole mergers. 
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there is no constant with the dimension of energy, it follows that the restriction of H^ DM to 
the static sector must vanish. The situation is quite different in Einstein- Yang-Mills theory 
where the Yang-Mills coupling constant g provides a scale. In c = 1 units, (gVG)^ 1 ~ Mass. 
Therefore, we can only conclude 

for some n-dependent constant C^ n \ As the horizon radius shrinks to zero, the static solution 
27| , |28[] under consideration tends to the solitonic solution with the same 'quantum numbers' 
n. Hence, by taking this limit, we conclude (gy/G)~ 1 C^ = M^ I ^ ) 1 {^ n ' ( • n ' , . Therefore, we have 
the following interesting relation between the black hole and solitonic solutions: 

where the integral of (3ft is evaluated on the 1-dimensional 'parameter space' of Sft (given by 
the horizon radius). Furthermore, as is clear from the above discussion, the ADM mass of 
the soliton is a multiple of {g\fG)~ 1 . Thus, somewhat surprisingly, the derivation of the first 
law in the isolated horizon framework has led to an interesting relation between the ADM 
masses of black holes and their solitonic analogs in the static sector. 



VIII. DISCUSSION 

In the first part of this paper, we introduced the notions of weakly isolated and isolated 
horizons. In contrast with earlier work ||[9,12|, the definitions allow for the possible presence 



of distortion and rotation at the horizon. In addition, the present definitions are more 
geometric and intrinsic; in particular, they never refer to a foliation. 

The notion of an isolated horizon, unlike that of an event horizon, is completely quasi- 
local. One can test if a given 3-surface in space-time is (weakly) isolated or not simply by 
examining space-time geometry at the surface. Furthermore, space-times admitting an iso- 
lated horizon A need not admit any Killing field even in a neighborhood of A. In particular, 
they can admit radiation in the exterior region. Therefore, such space-times can serve as 
more realistic models of late stages of a gravitational collapse or black hole merger. In ]TB| 



the near A geometry of vacuum solutions is examined in detail using similar techniques to 
those used at null infinity. The resulting structure — presence of a preferred 'rest frame', 
constrains on possible isometries, Bondi-type expansions of the metric — should be useful 
to extract physics in the strong field regime of general relativity, especially in the problem 
of binary black hole collisions. 

This paper, however, focused on another aspect of isolated horizons: extensions of the 
zeroth and first laws of black hole mechanics. All previous discussions of these laws were 
restricted to perturbations of stationary black holes. Using Lagrangian and Hamiltonian 
frameworks, we extended these laws to arbitrary space-times admitting non-rotating iso- 
lated horizons in Einstein-Maxwell-dilaton and Einstein- Yang-Mills theory. Furthermore, 
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the analysis suggests that it should be rather easy to incorporate other forms of matter, 
provided they admit Lagrangian and Hamiltonian descriptions. 

The generalization of black hole mechanics presented in this paper has several interesting 
features. First, all quantities that enter the first laws are defined locally at the horizon A. 
In standard treatments, some quantities such clS £1X6 cl and surface gravity are defined at the 
horizon. Others, like energy and sometimes [2v| even the Yang- Mills/Maxwell charge and 
potential, are evaluated at infinity. In part because of this 'mismatch', to our knowledge the 
'physical process version' of the first law [IIJ had not previously been established for processes 
which change the charge of the black hole. Since all quantities in the present treatment 
are defined locally at the horizon, it is now straightforward to establish the law for such 
processes ||. Secondly, other treatments based on a Hamiltonian framework [10,|32,23| often 
critically use the bifurcate 2-surface which does not exist in the extremal case. Therefore, 
extremal black holes are often excluded from the first law. The present analysis never makes 
reference to bifurcate surfaces (which do not exist in physical space-times resulting from 
gravitational collapse). Therefore, our discussion of the first law holds also in the extremal 
case. Thirdly, with obvious modifications of boundary conditions at infinity, our analysis 
includes cosmological horizons where thermodynamic considerations are also applicable [[II] . 



Finally, and perhaps most importantly, our analysis sheds new light on the 'origin' of the 
first law: it arose as a necessary and sufficient condition for the existence of a Hamiltonian 
generating time evolution. A new feature of our framework is the existence of an infinite 
family of first laws corresponding to the infinite family of 'permissible' vector fields t a . (A 
vector field t a is permissible if it is Hamiltonian, that is, induces canonical transformations on 
the phase space.) In theories where we have sufficient control on the space of static solutions, 
such as Einstein-Maxwell, one can select a natural evolution field t a Q . Corresponding to 
evolution along this t", there is a canonical notion of energy which can be interpreted as 
the mass of the isolated horizon. There exist also preferred values of surface gravity and 
electric potential and a canonical first law. This additional structure is extremely useful 
in other applications of the framework, such as extraction of physical information from 
numerical simulations of black-hole collisions. However, it is not essential to the discussion 
of mechanics: our derivation of the first law in Sections [V] and [VT] does not require any 
knowledge of the static sector of the theory. 

The Hamiltonian approach to black hole mechanics has appeared in the literature before, 
most notably in the work of Brown and York [[14]]. The spirit of the Brown- York approach 
is similar to ours. In particular, they do not restrict themselves to stationary situations. 
However, in that work, the focus is on an outer, time-like boundary whereas our focus 
is on the inner, null boundary representing the isolated horizon. Conserved quantities in 
presence of internal boundaries were recently discussed also by Julia and Silva pOl in a more 



general context of theories with gauge symmetries. As in our framework, their treatment 
exploits the simplifications that occur in a first order formalism and the final surface-integral 
expressions of conserved charges are dictated by the precise boundary conditions imposed at 
the internal boundaries. Their treatment is based on superpotentials and is thus complements 
the Hamiltonian methods used here and in [14]. 
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In this paper, the Lagrangian and Hamiltonian frameworks are based on real tetrads 
and Lorentz connections. It is therefore quite straightforward to extend our analysis to 
any space-time dimension. Indeed, it has already been extended to 2+1 dimensions in [ |3Tf . 
However, our phase space — and especially the explicit symplectic structure used here — is 
tailored to the Einstein-matter system. While it should be possible to extend it to higher 
derivative theories of gravity as in [02]| , that task would not be as simple. 
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APPENDIX A: EXAMPLES OF DISTORTED HORIZONS 



Because of the no hair theorems in the Einstein-Maxwell theory, distorted horizons have 
received a rather limited attention in the literature. Therefore, in this appendix we will 
discuss a few explicit examples in Einstein-Maxwell theory. For a general construction and 
an existence result, see |L8j,j33|]. 

To obtain explicit solutions, one has to impose symmetries. All solutions considered in 
this section will be static and axi-symmetric. As one would expect from the no-hair theorems, 
they fail to be asymptotically flat, whence they fail to represent isolated black holes in the 
standard sense. Nonetheless, they all satisfy the isolated horizon boundary conditions. That 
framework also serves to 'explain' the otherwise surprising feature that the surface gravity of 
these solutions depends only on the area and the charge and is insensitive to the distortion 
parameters. 

In the literature on static, distorted black holes, it is generally assumed that the solution is 
valid only in a finite region around the horizon and its distant behavior is suitably modified by 
the far-away matter which causes the distortion. For undistorted isolated horizons, Robinson- 
Trautman space-times [15] offer interesting examples of vacuum, asymptotically flat solutions 
which admit isolated horizons but no Killing fields whatsoever. Distorted analogs of these 
solutions are not known but presumably exist. It would be interesting to find them. 



1. Black Hole in a magnetic universe 



Let us begin with a simple example: an uncharged black hole in an 'external magnetic 
field' which distorts the horizon. The specific solution we wish to consider is static and 
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axisymmetric and was first obtained in the Ernst-potential framework ||34|| . The magnetic 
field is non-zero on the horizon. Thus, one has to consider the full set of Einstein-Maxwell 
equations on the horizon. 

The space-time has topology S 2 x 1R 2 and the metric is given by 



where 



2M\ , 2 dr 2 2 ,„„ 

dt 2 + — — + r 2 d6 2 

1 - 2M r 



r 2 sin 2 9 . , 2 . . „ , 

+ — f2~d<f> 2 (a.i; 



F = 1 + -B 2 r 2 sin 2 9. 
4 u 

Bq is a constant and on the axis the magnetic field is given by B = B$dr. Because F diverges 
at infinity, the metric fails to be asymptotically flat. For M = 0, the metric reduces to that 
of the Melvin universe and for M ^ it admits a Killing horizon at r = 2M. To examine the 
behavior at the horizon, let us first cast the metric in the Eddington-Finkelstein coordinates 
0, r, 9, <p) where dv = dt + (1 - dr: 

ds i = _F 2 fl--l dv 2 + 2F 2 dvdr + F 2 r 2 d9 2 + r ^^-d<p 2 . (A.2) 
V r J F l 

Since the metric is not asymptotically flat, the standard procedure of normalizing the Killing 
field to be unit at infinity is not applicable. Thus, we only have an equivalence class [i\ of 
(preferred) null normals to the horizon, £ oc J^. Let A be the Killing horizon and assume 
the Killing field ^ is a member of the equivalence class [£]. It follows trivially that (A, [£]) 
is a non-rotating isolated horizon. If Bq ^ 0, the scalar curvature 2 R of the horizon 2-metric 
has ^-dependence; the horizon is distorted. However, an explicit calculation shows that, as 
in the Schwarzschild space-time, the surface gravity k is given by l/2r = l/2i?A and the 
electrostatic potential $ vanishes on A. At first, it is quite surprising that while the presence 
of distortion affects /i, $n, ^2 and 2 R, it does not affect k or $. However, as we saw in 
Section 0, this result is to be expected from the general Hamiltonian considerations. 

2. Distorted black holes as special cases of Weyl solutions 

In this sub-section, we will review the construction of a large family of distorted black 



holes starting from Weyl solutions []35| and a recent generalization of these results to include 
electric charge |36fl . 

A general class of static, axisymmetric spacetimes was found by Weyl in 1917 P7|. The 
metric for such a spacetime can be cast in the following form: 

ds 2 = -e 2 ^dt 2 + e 2 ^ (dp 2 + dz 2 ) + e" 2 V^0 2 (A.3) 

where ip and 7 are smooth functions of p and z. Einstein's vacuum equations expressed in 
terms of ip and 7 take a particularly simple form. The equation for if}, 
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t PP + ~ + 4>,zz = 0, (A.4) 
P 

is simply the Laplace equation in flat space with cylindrical coordinates (p, z, </>). (In addition, 
i/j has to be independent of the angular coordinate (f>.) Given a solution for ip, the function 
7 can be determined by simple integration: 

l,z = 2/#>J ( A - 5 ) 

The Schwarzschild metric is of course a particular solution to these equations and corre- 
sponds to choosing for ip and 7: 

1, (L-M\ , 1, (L 2 -M 2 \ , A . 

W = V>s := 77 m r ■ and 7 = 73 := - In — 5- , (A.6) 



2 \L + M 

where L = |(/ + + 77 = - L) with /+ = yp 2 + (2 + M) 2 and L = ^p 2 + (z - M) 2 
and M is the mass of the Schwarzschild solution. Note that ip is just the Newtonian potential 
due to a rod of length 2M placed symmetrically about the origin on the z-axis. Both tps and 
75 diverge logarithmically in the limit p — > (for | z \ < M) . In order to recast this solution 
in the standard Schwarzschild form, one must transform from (z, p) to the Schwarzschild 
coordinates (r, 9) by 

z=(r-M) cos 9 and p 2 = r 2 (l - 2M/r) sin 2 (A. 7) 

This coordinate transformation shows that the horizon, r = 2M, corresponds to the line 
segment p = 0, | z \< M in Weyl coordinates. Therefore, the Weyl coordinates cover only 
the exterior of the horizon. 

Now, the key point is that ( |A.4j ), the only field equation one has to solve, is linear. Hence 



we can 'distort' the Schwarzschild solution simply by adding to an y solution ^ of the 
flat space Laplace equation which is regular along the z— axis [|35(| . Thus, we can set 

i) = i) S + 4>D and 7 = 7s + 7z> (A. 8) 



Substituting these expressions into ([A.4 ) and ( A.5f) and using the forms of the Schwarzschild 



functions ips an d 75, one can show that at p = 0, 

Id | p =o= 2^£> | p =o • (A. 9) 

This fact plays an important role in analyzing the horizon structure. 
In Schwarzschild coordinates, the distorted metric takes the form 

2(7 D -^ D ) 

rfs 2 = _ e 2 ^(l - 2M/r)dt 2 + dr 2 + e 2 ^ D ^ D) r 2 d9 2 + r 2 sin 2 0e~ 2fe # 2 . (A. 10) 

I J. £4 1\± 11) 
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As usual, the metric has a coordinate singularity at r = 2M. Let us therefore introduce the 
Eddington-Finkelstein coordinate v as before. The metric can be re-expressed in (v, r, 6, (ft) 
coordinates as 

d s 2 = -e 2 ^(l - 2M/r)dv 2 + (1 - 2M/r)- 1 e 2i>D (e 2 ^ D ~ 2 ^ - l)dr 2 

+ 2e 2 ^ D dvdr + e 2 ^~^\ 2 dQ 2 + e - 2 ^r 2 sin 2 Odcf 2 . (A.ll) 

Using condition (|A.9|) it is not difficult to show that the coefficient of dr 2 in the metric is 
regular at r = 2M ||35|| . 

It is immediately obvious from (|A.11|) that the r = 2M surface is a Killing horizon of Qr . 
However, we can not select a preferred normalization for this vector field since the metric is 
not asymptotically flat. As in the last sub-section, let A be the Killing horizon and choose 
£ oc J^. Then, it is straightforward to verify that (A, [£)) is a complete, non-rotating isolated 
horizon. Let us calculate the value of surface gravity for £ = We obtain 

k = (e 2 ^-^)— = — (A. 12) 

v J 2r 2r K ' 



a 



where we arrived at the last expression by using ( A.9|) . Again, while spin coefficient Re [fi] 



the Weyl component ^2 and the scalar curvature R of the horizon metric all depend on the 
distortion function ipD, somewhat surprisingly the surface gravity Km does not. 

The natural question is whether the above framework can be extended to obtain distorted 
Reissner-Nordstrom solutions. This turns out to be non-trivial because the key equation 
( |A.4j ) now acquires a source term from the electromagnetic field and this field itself depends 
non-trivially on ip through the Maxwell equations. At first, the coupled system appears to 
be hopelessly difficult. However, there exists a prescription |38[ for defining a new potential 
ifj in terms of ip and the electromagnetic field such that ip satisfies the flat space Laplacian 
( |A.4| ). Using this method, the known distorted black hole solutions were recently generalized 
to the charged case The distorted Reissner-Nordstrom solution is given by the metric 



ds 2 = -(1 - 2M/r + Q 2 /r 2 )e 2lpD dt 2 + (1 - 2M/r + Q 2 /r 2 )~ l e 2{ ~< D -^ D) dr 2 

+e 2{ lD -i, D ) r 2 d Q2 + e -2fe r 2 gin 2 M( p2 _ ^ ^ 

The forms of ip£, and 7^ are now substantially more complicated than in the uncharged case. 
Nonetheless, it is still possible to show that (|A.9|) continues to hold. As before this equality 



implies that the apparent singularity at r H — 2Mr H + Q = is only a coordinate singularity. 
The surface defined by r = rn is a Killing horizon of J^. There is once again, no natural 
way to normalize the Killing field, so we only have an equivalence class [£ a ] of null normals 
to the Killing horizon. (A, [£]) is a non- rotating isolated horizon. 
The surface gravity of 4 is given by 




(A.14) 
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Again, the surface gravity is independent of the distortion of the horizon and has the same 
dependence on the horizon radius i?A (which turns out to be equal to r H ) and charge Q as in 
Reissner-Nordstrom spacetime. Considerations of Section [V] suggest this peculiar behavior 
of k in all these examples is not accidental but can be 'explained' from general Hamiltonian 
considerations which led us to the first law. 



APPENDIX B: THE NEWMAN-PENROSE FORMALISM 
1. Notation and Conventions 



Let us begin with a summary of the Newman- Penrose formalism (see [^9| or |4CH42| for 



a complete account). Apart from the spacetime signature which we take to be (— , +, +, +), 
we will follow the conventions used in [Q. Consider a tetrad of null vectors n, £, m and m 
(n and £ are real while m is complex) which satisfy 

n.t = —1 n.m = n.m = , , 

i.m = £.rn = m.m = 1 . 

The directional derivatives along the basis vectors are denoted by 

D = tV a A = n a V a 5 = m a V a 5 = rrfV a . (B.2) 

The full the information contained in the connection is expressed in terms of twelve complex 
scalars called the Newman-Penrose spin coefficients defined as follows: 



(B.3) 



It is sometimes more useful to express these definitions in terms of covariant derivatives of 



K = 


-m a i b S7 b l a 


e = 


\ 


[fn a £ b V b m a - 


- n a £ b V b £ a ) 


7T 


= fn a £ b V b n a 


a = 


—m a m b V 'b£ a 


(3 = 


\ 


[m a m b Vbiri a 


- n a m b V b Q 


/' 


= m a m b V b n a 


P = 


—m a m b Vb£ a 


a = 


\ 


[m a m b 'VbTn a 


- n a fn b V b £ a ) 


A 


= m a m b V b n a 


T = 


-m a n b V b £ a 


7 = 


\ 


f 7jyf> n b^ b m a ■ 


- n a n b V b £ a ) 


V 


= m a n b V b n a 



the basis vectors: 












D£ = 


(e + e)£ - 


Km - 


- Km 


Dn = 


= — (e + e)n + nm + Wm 


A£ = 


(7 + i)t - 


Tm 


— Tm 


An = 


= — (7 + 7)71 + vm + Vm 


5£ = 


(a + (3)£ - 


■ pm 


— am 


5n = 


- — (a + (3)n + \im + Am 


Dm = 


T£ — nn + 


(e~ 


e)m 


Am = 


= vt — rn + (7 — 7)771 


5m = 


X£- an + 


((3- 


■ a)m 


5m = 


= ~p,£ — pn + (a — (3)m . 



(B.4) 



The ten independent components of the Weyl tensor are expressed in terms of five complex 
scalars \I/cb ^l; ^2, ^3 and ^4. The ten components of the Ricci tensor are defined in terms 
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of four real and three complex scalars $ cb ^11, $22, A, $ 10 , $ 2 o an d $21 • These scalars are 
defined as follows: 



*0 


= C abcd tmH c m d 


$01 = 


\Rabt a m b 


$10 = 


\R 


*1 


= C abcd i a m b i c n d 


$02 = 


^R a bin a m b 


$20 = 


\R 


^2 


= C a b c d£ a rn b m c n d 


$21 = 


\R ab rn a n b 


$12 = 


\R 


^3 


= C abcd tn b m c n d 


$00 = 


2 Rab^- ^ 


$ n = 


\R 


^4 


— C a b c dfn a n b m c n d 


$22 = 


\R ab n a n b 


A = 


R 

24 ' 



ad a m b 
ab m a m b 



(B.5) 





- 5(f) 




-50i 
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-50i 
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The six components of the Electromagnetic-Field 2— form F ab can be defined in terms of 
three complex scalars: 

o = -fWT ab 1 = _I(fV-m a m b )F ab 2 = n a m b Y ah . (B.6) 

The eight real Maxwell equations dF = and d*F = can be written as a set of four complex 
equations: 

(tt - 2a)0 o + 2p0x - k0 2 (B.7) 

-A0o + 27T0! + (p - 2e)0 2 (B.8) 

(2 7 - p)0 o - 2r0x + (B.9) 

z/0o - 2/i0! + (2/3 - r)0 2 . (B.10) 



2. Boundary Conditions 

In this section we describe the isolated horizon boundary conditions in the Newman- 
Penrose formalism. We will restrict ourselves to Einstein-Maxwell theory with zero cosmo- 
logical constant. 

In a null-tetrad adapted to the null hypersurface A, take £ to be a null normal, m and 
m tangent to A and n transverse to A. Since £ is hypersurface orthogonal and null, it is 
geodesic. This implies thatQ k np = and Im [p] = 0. Thus 

D£ b := £ a V a £ b = (e + e)t . (B.ll) 

The surface gravity is therefore given by = e + e and the expansion of £ is 9^ = Re [p] . 

For a non-expanding horizon A, the conditions on £ imply p = and the Raychaudhuri 
equation then implies a = and $00 = \R a b£ a £ b — 0. Furthermore, from ( |U.16| ) (which is 



16 



We will denote the NP spin coefficient k by k^p to distinguish it from the surface gravity n^y 
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a consequence of the energy condition), it follows that <f) = 0. This leads to the following 
conditions on the Ricci tensor at the horizon 



$00 = $ 01 = $io = 

$02 =0 $20 = $n = -2G0101 



(B.12) 



The first Maxwell equation (p.7[ ) gives 

Dfa = which implies D$n = . (B.13) 
Also, as shown in ( [11. 8| ) 

^o = and ^i = . (B.14) 

The intrinsically defined one-form uj a defined in ( |II.4|) is given by 

u a = -K(i)W a + (a + P)m a + (a + (3)m a . (B.15) 

It is often convenient to choose the null tetrad such that dn = which implies 

/i = ~p and n = a + j3 . (B.16) 

In this case we get a foliation of A spanned by m and m. Furthermore, by an appropriate 
spin transformation, we can choose e to be real so that e = I and thus the foliation is Lie 
dragged along £: 

C e m a = (e-e)m a = 0. (B.17) 

The one-form u now becomes 

u a = -K(e)n a + -nm a + W m a . (B.18) 

Let us consider a weakly isolated horizon (A, [£]). The condition Ctto = is equivalent to 
requiring 

£e7i~ = and C £ Km — (B.19) 

and as we proved in Section [LIB , these conditions imply that the surface gravity K(j) is 
constant on A. 

As mentioned in Section [II B| , a weakly isolated horizon with non-zero surface gravity 
admits a natural foliation. In the Newman- Penrose framework this foliation can be charac- 
terized as follows: It is the unique foliation on each leaf of which the pull-back of the 1-form 
7im a + 7f m a is divergence-free. This condition was first introduced by Hajicek J4!| in the 
context of stationary spacetimes. 

Finally, since our boundary conditions require that LiKm = 0, in a sense, a part of the 
zeroth law is simply assumed. As mentioned in Section |II Bj , we could have used a slightly 
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different set of boundary conditions which make no direct requirement on Km and yet lead 
to the zeroth law (as well as the results of Sections |T| - [VTJ). 



Let (A, [£]) be a non-expanding horizon, equipped with an equivalence class [£} of null 
normals to related to each other by constant positive rescalings. As above introduce a null 
tetrad where i is an element of [£} , m and m are tangent to the foliation, n is curl free and 
e is real. In place of Definition 2, let us assume that A admits a foliation by a family 
of 2— spheres transverse to [£} such that the NP spin coefficients in an associated null-tetrad 
satisfy: 

C e p = and £ e 7r = 0. (B.20) 

These conditions now replace the requirement C £UJ = used in the definition of a weakly 
isolated horizon. We can prove the zeroth law from these conditions as follows. First, 
consider the definition 2V[ a V&]£ c = R a bc d £,d of the Riemann tensor. In the NP formalism, 
these are written as a set of 18 complex equations known as the 'field equations'. For our 
purposes, we need only three of these equations W2 



Da — 5e = (p + 1 — 2e)a + (3a — (3e — k\ — K^y + (e + p)ir + $i 
D/3 - Se = (a + ir)a + (p - e)f3 - (p + 7)^ - (a - W)e + *i 
Dp - 5rr = (p - e - e)p + a\ + (W - a + 0)ir - uk + ^ 2 + 2A . 

Adding the first equation to the complex conjugate of the second equation and imposing our 
boundary conditions gives 

5( e + e) = 5n (e) = (B.21) 
while the third equation reduces to 

^ 2 = (e + e)p. (B.22) 

Equation ( B.21|) tells us that surface gravity is constant on each leaf of the foliation. It now 
only remains to show that it is also constant along I. To show this we turn to the Bianchi 
identity: Vr a i4c]<fe = 0- 111 t ne NP formalism, this is written as a set of nine complex and 
two real equations. We shall need only two of these equations pEfl 



Dm 2 - + A$ 00 - 5$oi + 2£>A 

= -A^ + 2(7r - a)*i + 3p^ 2 - 2/t^ 3 + a$ 02 

+ (2 7 + 27 - 7J)$ 00 - 2(a + r)$ 01 - 2r$ 10 + 2p$ u 

(B.23) 

D$n - 5$i + A$ 00 - 5$ i + 3DA 

= (2 7 + 27 - p - 7J)$ 00 + (tt - 2a - 2r)$ i + (W - 2a - 2r)$ 10 

+ 2(p + p)$n + a$ 02 + CT$20 - ^$12 - K$21 • 
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Subtracting these equations, imposing our boundary conditions and using A = 0, we get 
D^>2 — 0. Combining this result with (P.22j ) gives D(e + e) = 0. This completes the proof 
of the zeroth law within the alternate definition of weak isolation. Most of the results of this 
paper were first obtained using that definition. However, since that notion is tied so heavily 
to the presence of a foliation, its intrinsic meaning is somewhat obscure. Therefore, it was 
then replaced by Definition 2 used in the main body of the paper. 
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